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I  RESEARCH  WORK 


During  the  period  1  December  1981  -  31  January  1983  we  have  completed  a 
study  of  a  mechanism  which  can  cause  the  disruption  of  a  stationary  electric 
potential  structure.  Our  results  were  applied  to  the  parallel  electric  fields 
which  are  associated  with  auroral  arcs.  These  findings  were  reported  in  references 
1  to  A  listed  in  section  II.  Preprints  are  included  in  section  III.  This  work 
was  performed  in  collaboration  with  several  members  of  the  Department  of  Plasma 
Physics,  Royal  Institute  of  Technology,  Stockholm-Sweden. 

We  have  also  completed  a  program  of  research  directed  towards  analyzing  the 
structure  of  small  scale  auroral  vortices.  These  features  were  observed  by  the 
Air  Force  S3-2  satellite.  The  work  was  performed  in  collaboration  with  scientists 
at  AFGL.  Our  results  in  this  area  are  contained  in  references  5  and  7. 

At  this  point  we  shall  briefly  indicate  the  directions  of  future  research 
work  which  will  build  on  the  progress  made  during  the  past  five  years  of  this 
grant.  First  it  is  Intended  to  further  refine  our  model  describing  the  temporal 
development  of  a  plasma  sheath.  This  work  will  focus  on  the  mechanism  by  which 
the  dynamic  Bohm  criterion  is  established  at  the  sheath  edge.  Prof.  C.  Birdsall 
and  his  group  at  U.  C.  Berkely,  have  expressed  an  Interest  in  our  results  and  we 
anticipate  that  a  fruitful  collaborative  effort  will  be  established  in  this  area. 
Finally  we  will  continue  our  work  on  the  dynamic  properties  of  auroral  double 
layers  and  vortices.  It  is  hoped  that  our  previous  results  will  shed  some  light 
on  the  more  global  processes  that  are  associated  with  the  breakup  phase  of  a 
substorm.  Our  rationale  for  this  hypothesis  is  the  fact  that  the  small  scale 
features  that  we  have  focused  upon  are  signatures  of  this  epoch  of  the  storm. 

In  this  regard  we  have  submitted  a  proposal  to  AFGL  to  investigate  this  problem. 


At  present  the  study  has  been  approved  and  funding  should  commence  In  February 
1983. 

Before  concluding  this  section  we  would  like  to  thank  Dr.  H.  Radoski  for 
his  encouragement  and  support  of  our  research  efforts  over  the  period  of  this 


grant. 
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OBSERVATION  OF  SMALL-SCALE  AURORAL  VORTICES  BY  THE  S3-2  SATELLITE 

William  J.  Burke,  1 
Michael  Sllevltch,  2 
David  A.  Har<j|y,  * 

ABSTRACT 

We  have  studied  two  Intense  auroral  events  which  were  encountered  by  the  S3-2 
satellite  at  ~  0550  and  ~  1930  MLT,  during  a  substorm,  near  the  equatorward 
edge  of  the  region  1  current  system.  The  events  are  marked  by  large  deflections 
In  the  east -west  magnetic  field  component.  In  the  dawn  (dusk)  event  the 
deflection  was  100  (300)  nT  with  a  total  duration  of  4(7)  seconds.  In  both 
cases  the  deflection  corresponding  to  an  upward  current  sheet  was  sharpest, 
occurlng  over  .25  seconds.  Upward  current  strengths  of  45  and  135  u  A/m2 
within  latitude  extents  of  <  2  kn  are  inferred.  The  principal  region  of 
return  current  Is  not  contiguous  with  the  upward  current  sheet  and  is  more 
spatially  extended,  with  maximum  Intensities  In  the  range  of  10  to  15  u  A/m2. 

The  events  are  associated  with  electric  fields  whose  magnitude  can  exceed 
200  mV/m.  In  the  region  between  the  current  sheets,  adjacent  to  the  principal 
upward  current  sheet,  the  electric  field  rotated  by  180°  then  returned  to 
Its  original  orientation.  Such  electric  field  variations  give  rise  to  plasma 
vortex  flow  patterns  similar  to  these  observed  in  auroral  folds  and  curls. 
Oesplte  spatial  aliasing,  the  measurements  of  an  electron  detector  In  the 
vicinity  of  the  upward  current  sheets  provide  useful  information  concerning 
the  density  and  temperature  of  the  parent  populations  and  the  field-aligned 
potential  drop.  With  some  reasonable  assumptions.  It  Is  shown  that  observed 
values  of  j|  and  $|  are  consistent  with  collisionless,  single-particle 
theories. 

1.  Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA  01731 

2.  Electrical  Engineering  Department,  Northeastern  U.,  Boston,  MA  02115 


Introduction 

The  purpose  of  this  paper  Is  to  present  a  detailed  examination  of  two 
small-scale  electrodynamic  structures  observed  at  low  altitudes  by  the  S3-2 
satellite.  Both  events  were  detected  near  the  dawn-dusk  meridian.  In  the 
southern  auroral  Ionosphere,  during  a  substorm  period.  Large-scale  features 
of  the  substorm  period  have  been  simulated  by  the  Rice  University  group 
(Harel  et  al.,  1981  a,b;  Spiro  et  al.,  1981;  Karty  et  al.,  1982;  Chen  et  al., 
1982).  A  preliminary  analysis  of  electric  field  and  field-aligned  current 
measurements  during  the  events  has  been  given  by  Burke  (1981).  It  *  s  shown 
that  In  both  instances  the  satellite  passed  through  paired  current  »ets. 

The  upward  current  sheets  had  latitudinal  thicknesses  of  <  2  km  and  d  average 
current  densities,  j,  ,  of  45  and  135  u  A/m2.  The  return  (downwart  rent 
sheets  were  latitudlnally  more  extended  with  j,  In  the  10  to  15  y  A/m2  range. 
Also  the  principal  return  current  sheets  were  spatially  separated  from  the 
Intense  upward  current.  In  the  region  between  the  sheets,  but  close  to  the 
Intense  upward  currents,  the  electric  field  underwent  rotations  then  returned 
to  Its  original  orientation.  The  present  analysis  considers:  (1)1x1 
plasma  drift  motions  In  the  regions  of  electric  field  rotations,  and  (2) 
measurements  of  an  electron  spectrometer  In  the  vicinity  of  the  current 
sheets. 

Consideration  of  the  T  x  T  plasma  drift  motions  reveals  the  existence  of 
small-scale  vortices  embedded  In  regions  of  large-scale  sunward  plasma  drifts. 
The  vortices  are  qualitatively  similar  to  auroral  curls  observed  with  ground 
based  Instrumentation  (Halllnan  and  Davis,  (1970).  Based  on  concepts  first 
developed  by  Al f yen  (1950)  and  Webster  (1957),  Hall i nan  and  Davis  (1970) 
suggested  a  theoretical  model  for  auroral  curl  formation  where  vortices  result 
from  shears  In  the  ?  x  ff  drift  pattern.  This  feature  Is  Inherent  to  a  partially 
neutralized  auroral  charge  sheet.  Webster  and  Halllnan  (1973)  showed  that 


such  structures  are  unstable  so  that  Infinitesimal  ripples  distorting  the 
drift  velocity  grow  into  vortices.  The  process  Is  analagous  to  the  Kel vin-Helmholtz 
(K-H)  Instability  of  fluid  mechanics.  Recently  Wagner  et  al.  (1981)  have 
demonstrated  by  a  computer  simulation,  how  the  electrostatic  K-H  mechanism 
leads  to  the  formation  of  vortices  in  the  auroral  charge  sheets. 

Two  essentially  different  theoretical  approaches  to  the  relationship 
between  field-aligned  currents  (j j)  and  potential  ($|).  The  first 
approach  concerns  the  transmission  of  electrical  Information  between  magnetospheric 
generators  and  ionospheric  loads  via  kinetic  Alfven  waves  (Goertz  and  Boswell, 

1  979;  Lysak  and  Carlson,  1981;  Lysak  and  Dum,  1982).  The  second  approach 
assumes  a  potential  distribution  along  magnetic  field  lines  and  calculates  the 
field-aligned  current  reaching  the  ionosphere  (Knight,  1973;  Lemaire  and 
Scherer,  1973,  1974;  Whipple,  1977;  Fridman  and  Lemaire,  1980;  Lyons ,  1981). 

The  relationship  was  derived  by  considering  collisionless,  single-particle 
trajectories  in  a  magnetic  field.  Values  of  at  ?  « nhanced  by  $g  which 
opens  the  atmospheric  loss  cone.  The  degree  to  which  the  loss  cone  is  opened 
depends  on  the  altitude  and  the  magnitude  of  *j.  The  derived  relationship 
between  j  g  and  has  been  shown  to  be  consistent  with  rocket  measurements 
of  precipitating  electron  fluxes  above  an  auroral  arc  ( Lyons ,  1981).  In 
this  paper  we  compare  electron  flux  measurements  with  predictions  of  the 
second  model.  Validation  of  the  first  model  is  beyond  the  capability  of 
S3-2  Instrumentation. 

The  following  section  provides  a  brief  description  of  S3-2  instrumentation 
and  a  detailed  presentation  of  measurements  in  the  vicinity  of  the  two  auroral 
events.  For  the  sake  of  coherence,  some  aspects  of  the  measurements  previously 
reported  by  Burke  (1981)  are  repeated  here.  The  discussion  section  focuses 
on  how  the  observations  compare  with  the  theories  of  auroral  curls  and  j  |/*f 
relationships. 


Instrumentation 


The  S3-2  satellite  was  launched  Into  polar  orbit  during  December  1975 
with  an  initial  apogee,  perigee  and  inclination  of  1557  km,  240  km  and  96,3°, 
respectively.  It  is  spin-stabilized  with  a  nominal  spin  period  of  20  sec. 

The  spin  axis  is  nearly  perpendicular  to  the  orbital  plane.  The  orbital 
plane  itself  is  subject  to  a  slow  westward  precession  of  1/4°  per  day.  The 
scientific  package  on  S3-2  includes:  (1)  an  electric  field  experiment  consisting 
of  two  dipoles,  (2)  a  trlaxial  fluxgate  magnetometer,  and  (3)  an  energetic 
electron  spectrometer.  Electric  and  magnetic  field  components  are  usually 
calculated  in  a  spacecraft  centered  coordinate  system  with  X  positive  along 
the  satellite  velocity,  2  is  positive  toward  local  nadir  and  Y  completes  the 
right  hand  system.  When  the  orbital  plane  is  near  the  dawn-dusk  meridian, 
as  in  the  cases  studied  here,  Y  is  positive  in  the  antisunward  direction 
(cf.  Figure  1  of  Smiddy  et  al. ,  1980).  In  this  paper,  data  are  also  presented 
in  geomagnetic  coordinates. 

The  instrumentation  of  S3-2  and  the  methods  used  to  reduce  data  have 
been  described  in  detail  by  Burke  et  al.  (1980).  Briefly,  one  of  the  dipoles 
has  a  length  of  29.72  m.  It  lies  in  the  satellite’s  spin  plane  and  measures 
the  electric  field  component  along  the  satellite  trajectory  (Ex).  The  second 
dipole  has  a  length  of  11.18  m  and  lies  close  to  the  satellite's  spin  axis. 

As  discussed  by  Hozer  et  al .  (1980),  the  spin-axis  aligned  dipole  frequently 
experiences  dc  offset  problems  so  that  only  the  spin  plane  data  are  normally 
presented.  The  problems  arise  from  variations  in  the  plasma  density  and/or 
chemical  composition.  The  variations  in  dc  offset  due  to  changes  in  plasma 
characteristics  usually  occur  on  time  scales  much  longer  than  <  10  sec  required 
for  the  satellite  to  traverse  the  intense  electrodynamic  shock  structures 
studied  below.  We  have  examined  simultaneous  measurements  from  two  driftmeter 
experiments  on  S3-2  (8urke  et  al.,  1980)  and  have  found  no  evidence  of  small-scale 
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variations  In  the  Ionospheric  plasma  characteristics  while  crossing  the  structures 


In  such  Instances  the  spin-axis  component  of  the  electric  field  (Ey)  can  be 
determined  by  subtracting  dc  constar  offsets.  Potential  differences  between 
the  ends  of  the  dipole  are  sampled  22  times  per  second.  To  resolve  spatio-temporal 
structures  with  the  highest  possible  resolution,  all  of  the  electric  field 
measurements  are  used.  Magnetic  field  measurements  were  taken  at  a  rate  of  32 
samples  per  second,  with  a  one  bit  resolution  of  5  nT.  All  of  the  data  are 
presented  as  differences  between  the  measured  and  IGRF  75  model  components. 

The  electron  spectrometer  measures  electrons  with  energies  between  80  eV  and 
17  KeV  In  32  energy  channels,  compiling  a  full  spectrum  every  second.  The 
aperture  of  the  instrument  is  mounted  in  the  spin  plane  of  the  satellite. 

Clearly,  the  particle-characteristics  of  structures  through  which  the  satellite 
passes  In  0.25  sec  cannot  be  resolved  fully  with  this  detector.  As  seen  below, 
through  happy  coincidences  of  energy-sampling  and  look  direction,  some  useful 
information  about  precipitating  electrons  has  been  retrieved. 

EXPECTED  SIGNATURES  OF  AURORAL  VORTICES 


Since  this  paper  largely  concerns  auroral  vortices  it  Is  useful  (before 
examining  the  data)  to  consider  their  expected  signatures  In  S3-2  measurements. 
Here  and  throughout  this  paper  we  assume  that  vortices  are  time-stationary  over 
the  Interval  of  satellite  encounter  and  are  reflected  in  the  T  x  F  drift  motion 
of  Ionospheric  plasma.  This  differs  somewhat  from  auroral  vortices  studied  by 
optical  means.  Optical  signatures  reflect  the  characteristics  of  energetic, 
auroral  electrons.  Such  electrons  may  have  more  complex  drifts  than  that  due 
to  I  x  I  motion.  Also,  If  there  are  field-aligned  potential  drops  the  F  x  F 
drift  of  energetic  electrons  and  conseguent  optical  vortlclty  may  be  larger  than 
the  vortlclty  of  cold  plasma  at  S3-2  altitudes. 

With  these  cautions  In  mind  let  us  consider  the  situation  sketched  in 
Figure  1,  which  Is  designed  to  be  representative  of  the  southern  hemisphere  In 
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the  early  mornlnq  sector.  For  simplicity  assume  that  the  satellite  moves 
southward  with  a  velocity  7S  along  a  magnetic  meridian.  The  magnetic 
field  Is  directed  radially  outward  from  the  earth  along  the  -Z  axis.  The 
satellite  trajectory  carries  It  to  the  west  of  a  flux  tube  containing  a 
small,  negative  space  charge.  As  shown  In  Figure  1  the  electric  field  due 
to  the  negative  charge  Tc  Is  directed  radially  toward  the  flux  tube.  The 
cold  plasma  has  a  counter  clockwise  rotation  7C  »  (Fc  x  1T)/B ^  around  the 
charged  flux  tube.  The  top  panel  of  Figure  1  shows  that  if  the  only  electric 
field  were  due  to  the  space  charge,  the  satellite  would  measure  a  field  with 
an  X  component  along  the  satellite  trajectory  while  the  satellite  was  to  the 
north  of  the  charge.  This  component  would  go  to  zero  at  the  point  of  closest 
approach  and  reverse  direction  to  the  south  of  the  charged  flux  tube.  The 
potential  obtained  by  integrating  Tc  along  the  satellite  trajectory  (second 
panel) 

•c  ■  -  ft c  -Vs  dt 

has  a  minimum  value  at  the  point  of  closest  approach.  The  potential  due  to  a 
postlvely  charged  flux  tube  would  have  a  maximum  at  the  point  of  closest  approach. 

The  morning  sector  of  the  auroral  oval  is  marked  by  an  equatorward- 
directed,  large-scale  convection  electric  field  (TT0)  that  has  a  relatively 
small  westward  component.  The  T0  x  IT  drift  (70 )  is  mostly  In  the  sunward 
(eastward)  direction.  The  potential  obtained  by  Integrating  E0  along  the 
satellite  trajectory  (third  panel  of  Figure  1) 

«b  *  -  /  f0  •  7S  dt 

Is  a  monotoically  Increasing  function. 

The  trace  In  the  bottom  panel  of  Figure  1  gives  the  total  potential  *s 
measured  by  the  satellite,  the  superposition  of  *c  *o»  It  has  a  local 
maximum  to  the  north  of  the  space  charge  and  a  local  mlnlnum  near  the  point 
of  closest  approach.  Since  cold  plasma  drifts  along  equlpotentlal  contours. 


3 


streamlines  of  the  flow  are  obtained  by  connecting  equl potential  points. 

The  direction  of  equipotential  contours  are  obtained  from  the  total  velocity 
T0  +  VC  of  the  plasma.  Mite  that  In  the  case  presented  In  Figure  1  a  vortex 
In  the  rest  frame  of  the  plasma  appears  as  an  S-like  reversal  of  the  plasma 
flow  In  the  satellite  frame  of  reference. 

OBSERVATIONS 

He  now  consider  the  morphology  of  the  two  electrodynamic  structures. 

They  were  detected  on  opposite  sides  of  the  oval  during  a  substorm  period  on 
19  September  1976.  In  analyzing  data  we  have  assumed  that  S3-2  encountered 
time-stationary,  spatial  structures.  That  Is,  the  time  scales  for  changes  in 
the  quasi-dc  electric  and  magnetic  fields  is  long  in  comparison  with  the 
satellite  traversal  period.  Attention  was  first  drawn  to  these  events  by 
their  magnetic  rather  than  electric  field  signatures.  Figures  2a  and  2b  are 
plots  of  a  By  measured  on  the  dawn  and  dusk  side,  southern  hemisphere 
auroral  oval  passes  of  Rev  4079A  and  4079B,  respectively.  The  data  are 
plotted  as  functions  of  UT  (in  seconds  of  the  day  and  in  hours  minutes), 
invariant  latitude,  magnetic  local  time  and  altitude.  The  first  event  begins 
at  0953:06  UT  ( A  *  68°,  MLT  =  05.4),  and  appears  as  a  100  nT  excursion  of  A  By 
that  lasts  for  4  seconds.  The  second  event  is  a  300  nT  excursion  of  A  By 
beginning  at  1148:16  UT  (A  *  67°,  MLT  a  19.3)  and  lasting  for  7  seconds. 

The  excursions  are  from  a  baseline  established  by  the  large-scale  region  1 
and  region  2  current  system.  In  both  cases  the  excursions  are  located  near 
the  equatorward  edge  of  the  region  1  current  sheet.  The  altitudes  of  observation 
were  1050  and  1350  km,  respectively.  The  Blrkeland  current  densities  associated 
with  these  A  By  excursions  greatly  exceed  those  normally  found  near  auroral 
arcs  (Choy  et  al.,  1971). 

A  detailed  analysis  of  the  events  was  performed  by  examining  the  individual 
electric  and  magnetic  field  measurements  across  the  events.  Rapid  fluctuations 
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(not  shown)  are  found  In  the  potential  differences  measured  across  both 
spinning  and  axial  dipoles.  The  fluctuations  are  Qualitatively  similar  to 
those  measured  at  higher  altitudes  near  small-scale  electrostatic  "shocks" 

(Mozer  et  al . .  1977).  The  satellite  was  sufficiently  high  In  altitude  that 
variations  In  the  axial  dipole  measurements  across  the  events  were  attributed 
to  changing  Ey  rather  than  to  changing  DC  offsets  on  the  active  elements  of 
the  electric  field  sensors.  Values  of  Ey  were  calculated  by  subtracting 
the  electric  field  component  due  to  satellite  motion  across  magnetic  field 
lines  and  our  best  estimate  of  the  DC  offset.  There  may  still  remain  a 
small,  unknown  DC  bias  In  the  measurements,  however,  relative  variations 
in  Ey  are  accurate.  The  electric  field  components  were  transformed  from 
satellite  to  geomagnetic  north-south  (Efg)  and  east-west  (E^y)  coordinates 
by  assuming  that  T  •  IT  »  0  everywhere.  The  three  components  of  aE  were 
also  computed  in  geomagnetic  coordinates.  For  the  remainder  of  this  section 
we  use  the  highest  possible  time  resolution  for  S3-2  measurements  of  T  and  "3T, 
considering  first  the  dawn-side  then  the  dusk-side  structure. 

(i )  Dawnside  Structure 

Figure  3  is  a  plot  of  four  and  one  half  seconds  of  electric  field  measurements 
In  geomagnetic  coordinates  starting  at  35586.2  (0953:06.2)  UT.  Current  densities, 
derived  from  the  slopes  of  A  Bew,  are  represented  in  histogram  form  with  positive 
values  corresponding  to  currents  out  of  the  ionosphere.  A  Byo  and  a  Bf5  show 
no  systematic  variations  across  the  region  of  Interest.  In  the  quarter  of  a 
second  after  35586.9  UT  a  Bew  increased  from  350  to  450  nT.  There  is  a 
partial  recovery  In  the  next  eighth  of  a  second  followed  by  a  1.25  second 
period  where  A  Bew  shows  a  very  weak  slope.  This  Is  followed  by  a  pair  of  current 
sheets.  The  principal  recovery  occurs  In  approximately  3/4  of  a  second  after 
this  pair  of  current  sheets  .  The  most  Intense  current  density  Is  45  y  A/m2 
directed  out  of  the  Ionosphere.  The  return  current  is  more  dispersed  over 
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the  event.  The  latitudinal  thickness  of  the  upward  current  sheet  projected 
to  an  altitude  of  100  km  Is  ~  1.5  km,  a  typical  width  of  optical  arcs  when 
viewed  from  the  ground  In  the  direction  of  magnetic  zenith. 

Vector  representations  In  the  bottom  panel  show  that  the  electric  field 
was  originally  oriented  to  the  west  of  north  and  remained  that  way  Into  the 
Intense  current  sheet.  Between  35586.2  and  35587.0  UT  Efy  had  a  nearly 
constant  westward  value  of  ~  20  mV/m.  During  the  same  interval  the  northward 
component  Increased  from  50  to  200  nV/m.  In  the  dawnslde,  southern  hemisphere 
an  Increasing  northward  electrostatic  field  component  In  satellite  coordinates 
corresponds  to  a  region  of  negative  electric  field  divergence.  Thus,  this 
upward  current  sheet  is  also  a  region  of  negative  space  charge.  After  35587.2 
the  electric  field  rotated  through  east  to  southeast  then  back  through  east 
to  northwest.  It  is  Important  to  note  that  Efy  and  Ej$  underwent  major 
reversals  in  a  region  where  "aTTwas  essentially  constant. 

Figure  4  Is  an  expanded  plot  of  &  B^y  and  E^  (top  two  traces)  In  the  vicinity 
of  the  upward  current  sheet.  The  position  of  the  upward  current  sheet  is  represented 
by  the  shaded  area  In  the  bottom  panel.  The  potential  $  (third  trace) 

Is  obtained  by  Integrating  the  electric  field  along  the  satellite  trajectory. 

The  potential  has  two  extrema  with  *  =  0  arbitrarily  set  at  the  poleward 
extremum.  Every  fifth  electric  field  vector  In  Figure  3  and  the  equi potential 
points  marked  X  and  0  has  been  used  to  illustrate  the  plasma  flow  velocity  V  In 
the  fourth  panel  of  Figure  4.  Since  the  T  x  ¥  plasma  flow  Is  directed  along 
equi potential  surfaces  plots  of  *  and  V  can  be  combined  to  approximate  the  two 
dlmentlonal  flow  characteristics  of  plasma  In  the  vicinity  of  the  arc  (bottom 
panel).  In  drawing  the  flow  lines  we  have  endeavored  to  maintain  a  divergence 
free  flow  pattern  In  two  dimensions  since  V  •  (T  x  ¥)  vanishes  in  the  steady 
state  (If  T  .¥  -0).  Comparing  the  potential  distribution  In  Figure  4  with 
the  bottom  panel  of  Figure  1  we  note  that  the  potential  maximum  at  0953:06.2  UT 
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lies  well  above  the  projection  of  the  potential  due  to  the  convection  electric 
field.  This  suggests  that  the  satellite  passed  to  the  east  of  a  postlvely  charged 
flux  tube.  The  point  of  closest  occurred  at  0953:06.2  UT.  The  satellite  than 
pased  to  the  west  of  a  negatively  charged  flux  tube  with  closest  approach  occurring 
at  0953:06.6  UT.  When  viewed  antiparallel  to  ¥  In  the  rest  frame  of  the  plasma 
one  finds  two  vortices  with  clockwise  and  counterclockwise  rotations  around  the 
positively  and  negatively  charged  flux  tubes,  respectively.  Assuming  a  radius 
of  ~2  km  for  the  counterclockwise  vortex  we  estimate  the  vortlclty  to  be 
~1.2  sec-^.  As  shown  in  Figure  4,  the  counter-clockwise  vortex  Is  detached 
from  both  the  upward  auroral  current  sheet  and  the  adjacent,  weak  return  current 
sheet. 

Since  the  S3-2  electron  spectrometer  requires  a  full  second  to  compile  a 
32  point  energy  spectrum  care  must  be  exercised  interpreting  measurements 
in  the  vicinity  of  the  spatially  narrow  current  sheet.  Figure  5  provides  five 
sequential  measurements  of  directional,  differential  flux.  Dashed  lines  Indicate 
the  one  count  per  sampling  interval  flux  level.  To  the  right  of  the  figure  are; 

(1)  the  UT  at  which  the  spectrum  compilation  began  and,  (2)  the  pitch  angle 
(  c$  sampled  half  way  through  the  second.  In  compiling  spectra,  energy  channels 
were  alternately  sampled  in  Increasing  then  decreasing  energy.  Energy  channels 
were  sampled  from  high  to  low  energies  in  the  cases  of  spectra  1,  3  and  5  and 
from  low  to  high  energies  in  the  cases  of  spectra  2  and  4. 

Spectra  1  and  2  show  broad  thermal  distributions  typical  of  diffuse  auroral 
precipitation.  Spectrum  3  Includes  electrons  with  pitch  angles  between  108°  and 
90°.  Like  spectra  1  and  2,  high  energy  channel  measurements  were  at  or  near 
the  1  count  level.  There  Is,  however,  an  electron  burst  with  a  peak  energy  of 
~3  keV  that  coincides  exactly  with  passage  through  the  Intense  upward  current 
sheet.  As  the  detector  turned  toward  nadir,  (spectra  4  and  5)  count  rates 
decreased  then  vanished.  After  considering  the  duskslde  structure,  we  return 
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to  discuss  the  electron  burst  of  spectrum  3  and  a  similar  duskslde  precipitation 
event  In  more  detail. 

(11 )  Duskslde  Structure 

The  electric  and  magnetic  field  components  plotted  In  Figure  6  show  that 
the  duskside  structure  had  the  more  complex  structure.  Near  the  poleward  boundary 
of  the  structure,  Efg  decreased  from  75  to  ~Q  mV/m  between  42496.8  and  42497.6  UT 
while  Efy  maintained  a  nearly  constant,  low  value.  After  this  time,  both 
components  Increased  In  magnitude  then  reversed  directions.  The  magnitude  of 
E  after  the  reversal  exceeded  300  mV /pi.  The  components  then  returned  to  their 
original  values.  After  42498.0  UT  Ef$  underwent  large  variations  in  magnitude, 
but  generally  remained  southerly  in  orientation.  E^y  underwent  large  varia¬ 
tions  in  both  magnitude  and  direction  before  returning  to  a  near  zero  value  at 
42504.8  UT. 

Magnetic  field  variations  across  the  evening  side  structure  show  patterns 
that  on  a  large  scale  are  similar  to  those  of  the  morning  side  structure.  In 
the  quarter  of  a  second  prior  to  42497.8  UT,  4  B^y  decreased  by  300  nT,  corresponding 
to  an  average  current  density  over  this  interval  of  135  u  A/m2.  As  in  the  previous 
case  there  is  a  weak  return  current  immediately  adjacent  to  the  upward  current  but 
the  principal  recovery  in  4  B^y  occurs  in  a  region  spatially  seperated  from  the 
upward  current.  Again,  the  electric  field  underwent  a  rotation,  this  time  of 
360°  at  the  equatorward  edge  of  the  intense  upward  current  sheet  and  weak  return 
current  sheet  pair.  The  Internal  structure  of  the  evening  side  event,  however. 

Is  much  more  complex  than  that  of  the  morning  side  event.  This  Is  reflected  In 
the  "multiple  current  sheets"  found  between  42497.8  and  49500.8  UT.  The  Intensities 
of  these  "current  sheets"  were  calculated  from  slope  of  4  B^y  i.e.  assuming  that 
the  satellite  crossed  a  longitudinally  extended  structure  at  normal  incidence. 

The  presence  of  significant  variations  In  4  Bfg  as  well  as  4  Bfy  Indicates 
that  the  "multiple  current  sheet"  level  of  complexity  Is  too  simple. 
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If  aB£u  and  a  B^g  varied  In  the  same  sense  across  the  structure,  then  the 
measurements  would  Indicate  that  the  satellite  trajectory  crossed  a  series  of 
current  sheets  at  an  angle  other  than  90°.  Indeed  similar  variations  of  the  two 
components  near  the  poleward  edge  of  the  structure  Indicate  that  the  arc  was 
crossed  at  an  angle  of  9°  away  from  normal  Incidence.  The  anti  correlated 
variations  In  aBj^  and  a  Bgy  near  42500  UT  suggest  that  $3-2  passed  close 
to  a  line  current  embedded  within  an  extended  sheet  current. 

A  different  complexity  Is  found  In  the  measurements  taken  between 
42498.8  and  42499.8  UT.  Variations  In  II  coincide  with  variations  in  the  mag¬ 
nitude,  but  not  direction,  of  F.  The  measurements  are  consistent  with  the  satellite 
having  detected  large  amplitude  Alfven  waves.  The  relationship  between  the  ampli¬ 
tudes  of  the  electric  and  magnetic  fields  of  Alfven  waves  comes  directly  from 
Maxwell's  equation;  6  E  =  V/^  6  B.  The  Alfven  speed,  V^  *  B0/  /p0p  ,  where  B0 

is  the  main  field  of  the  earth,  and  p  the  mass  of  the  plasma.  In  the  centered 

dipole  approximation,  at  an  altitude  of  1350  km  along  the  A  =  67°  field  line 
B0  *3.3  x  10"5  T.  For  a  dominant  ion  species  of  0+  and  a  number  density  between 
1010  and  109  m-3 ,  is  between  1.8  and  5.7  x  106  m/sec.  Thus,  for  an  average 

variation  in  6  B  of  75  nT,  as  seen  in  the  period  of  interest,  6  E  should 

vary  between  135  and  420  mV/m.  The  variation  in  F  during  this  time  is  ~  300 

mV/m.  Ten  seconds  prior  to  encountering  the  structure  a  planar  thermal  ion 
sensor  on  S3-2  measured  a  plasma  density  of  3  x  lO^tr3.  Note,  however, 
that  a  1  Inear-regime  Alfven  wave  can  not  explain  the  variations  in  A  B^y 
near  42497.8  UT.  A  300  nT  magnetic  field  amplitude  would  require  a  6  E 
between  540  and  1700  mV/m,  far  greater  than  the  actually  mesured  elec¬ 
tric  field.  In  the  non-linear  regime  however,  Lysak  and  Dum  (1982)  have  shown 
that  such  a  simple  relationship  between  aB  and  5E  no  longer  maintains. 

Finally,  we  note  that  even  If  the  variations  of  E  and  A  B  between  42498.8 
and  42499.8  UT  are  due  to  Alfven  waves,  It  Is  not  possible  to  specify  directions 
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of  propagation  uniquely.  Since  E  varied  In  magnitude,  but  not  direction, 

$7  lies  approximately  In  the  1C-SW  direction.  Unfortunately,  establishing 
an  "unperturbed  baseline"  from  which  we  measure  $7  demands  an  unacceptably 
high  level  of  subjectivity. 

Figure  7  is  an  expanded  plot  of  the  potential  (top  panel)  and  flow  vectors 
(bottom  panel)  In  the  Immediate  vicinity  of  the  electric  field  rotation  shown  in 
Figure  6.  As  In  the  case  shown  In  Figure  4  the  potential  has  two  extrema.  In 
this  case,  however,  the  potential  maximum  lies  close  to  the  average  convection 
potential.  Even  in  the  satellite  frame  of  reference  the  flow  is  dominated  by 
a  counterclockwise  rotation  about  a  negatively  charged  flux  tube.  The  point 
of  closest  approach  occurred  at  42497.8  UT.  We  estimate  the  radius  of  the  vortex 
to  be  0.75  km  and  *  2.5  sec-1. 

Figure  8  provides  six  consecutive  spectral  measurements  taken  in  the 
vicinity  of  the  duskslde  event.  Odd  (even)  numbered  spectra  were  sampled 
from  low  (high)  to  high  (low)  in  energies.  Poleward  of  the  event  (spectrum  1), 
the  electron  flux  was  close  to  or  below  instrumental  sensitivity.  This  is 
also  true  of  the  highest  energy  channels  sampled  in  spectrum  2.  Spectrum  2 
is  marked  by  a  burst  of  electrons  with  a  peak  flux  of  8  x  109e/(cm2- 
sec-ster-keV)  at  3.5  keV  and  a  *  151°.  The  slopes  of  the  low-energy 
portions  of  both  spectra  2  and  3  are  relatively  flat.  A  peak  flux  5  x  109  e/( cm^ 
sec-ster-keV)  at  1.7  keV  was  measured  In  spectrum  3  while  the  detector 
was  looking  close  to  the  direction  of  the  magnetic  field.  Following  the 
broad  thermal  distribution  of  spectrum  4,  the  flux  rapidly  decreased  to 
background.  Obviously  there  Is  a  high  degree  of  spatial -temporal  aliasing 
In  spectra  2  and  3.  However,  two  purely  phenomenological  remarks  can  be 
made.  First,  a  detailed  examination  of  magnetometer  and  particle  detector 
data  records  shows  that  the  3.5  keV  electron  burst  of  spectrun  2  was  encountered 
1/4  sec  prior  to  (2  km  poleward  of)  the  intense  upward  current  sheet.  Second, 


a  flux  of  1.2  x  10l°e/(cm2-sec-ster-keV)  for  electrons  with  energies  near  80  eV, 
seen  In  spectrum  2,  Is  the  highest  level  observed  for  $3-2 's  lowest  energy 
channel. 

It  Is  possible  to  glean  further,  albeit  somewhat  speculative.  Information 
by  plotting  spectra  3  and  2  of  Figures  5  and  8,  respectively,  as  distribution 
functions,  f  In  Figure  9.  The  distribution  functions  of  both  the  morning  and 
evening  side  electron  bursts  have  peaks  at  energies  of  3.0  and  3.5  keV. 

Following  Burch  et  al.  (1976)  we  Interpret  electrons  with  energies  above  the 
peak  as  accelerated  primaries.  Those  with  energies  below  the  peak  are  interpreted 
as  secondary  and  degraded,  primary  electrons.  In  our  semiloq  plot  of  f  versus 
energy,  the  points  at  and  beyond  the  peak  have  been  fit  to  straight  lines. 

The  burst  electrons  come  from  magnetospheric  populations  that  were  Maxwellian 
In  energy  distribution  and  accelerated  through  field-aligned  potential  drops 
of  3.0  and  3.5  keV.  The  slopes  of  the  straight-line  fits  indicate  magnetospheric 
temperature  of  ~  200  eV.  The  peak  value  of  the  distribution  function  f0  can 
be  used  to  estimate  the  densities  of  the  parent  populations.  In  convenient 
units, 

(1)  f0(cm~6-sec"3)  =  8.61xl0“2^  n(cm"3)/T*,5(ev). 

For  the  morning  side  f0  was  estimated  by  assuming  a  potential  drop  of  3.0  kV  and 
then  extrapolating  the  fit  to  3.0  keV.  This  gives  a  value  of  ~  9  x  10-29 
(cm_6-sec“3).  For  the  evening  side  fQ  was  estimated  by  determining  the 
intersection  of  the  linear  fit  above  the  peak  and  the  linear  fit  for  the 
three  points  below  the  peak.  This  gives  a  value  of  ~2.8  x  10"28  (cm-6 -sec”3). 
Assuming  magnetospheric  Isotropy  we  estimate  the  parent  densities  to  be 
0.30  cm“3  and  0.92  cm-3  for  the  morning  and  evening  side  events,  respectively. 
These  temperatures  Indicate  that  the  events  occur  In  the  boundary  plasma  sheet 
as  defined  by  Winnlngham  et  al.  (1976).  In  the  evening  side  event  there  is 


also  evidence  for  a  higher  energy  population  whose  temperature  is  4  keV  and 


density  Is  1.1  an~3.  For  simplicity  In  later  calculations  Me  assume  that 
this  population  has  also  been  accelerated  through  a  3.5  kV  potential 
drop.  Anticipating  The  discussion  section  we  compute  the  quantity 
(2)  j0  *  nq  /  kf/2jr  m 

for  each  of  these  populations.  This  quantity  Is  the  field-aligned  current  density 
that  would  be  produced  by  an  Isotropic  electron  distribution  in  the  absence  of 
a  field-aligned  potential  drop.  Table  1  gives  a  compilation  of  the  n,  kT,  *| 
and  j0  calculated  from  Figure  9. 


TABLE  1 


Event 

n{cm"3) 

kT(eV) 

(kV) 

jn[  I*/""2 

Morning  (Cold) 

0.30 

200 

3.0 

1.2 

Evening  (Cold) 

0.92 

200 

3.5 

3.45 

(Hot) 

1.1 

4000 

3.5 

18.6 

Discussion 

In  this  section  we  discuss  the  S3-2  measurements  In  the  light  of  our 
present  theoretical  understanding  of  (1)  the  generation  of  auroral  vortices 
and  (2)  the  relationship  between  field-aligned  currents  and  potential 
drops.  Although  other  aspects  of  the  event  are  of  importance,  we 
restrict  ourselves  In  this  paper  to  considering  only  these  two  points. 

Advances  In  the  technology  of  low-light-level  television  cameras  and  other 
space-optical  systems  have  allowed  the  classification  of  many  dynamical  auroral 
features  (Qqutl ,  1981  and  references  therein).  At  least  three  classes  of 
auroral  vortices  have  been  Identified:  spirals,  folds  and  curls.  Spirals  are 
large-scale  vortices  with  diameters  in  the  20  to  1300  km  range  (Davis  and 
Halllnan.  1976).  When  viewed  antlparallel  to  the  magnetic  field  j?  the 
flow  pattern  Is  clockwise.  Halllnan  (1976)  argued  that  spirals  result  from  a 
Kel vln-Helmholtz  like  current  sheet  Instability.  Both  folds  and  curls  are 


small-scale  phenomena  (Hallinan  and  Davis,  1970)  Folds  have  wavelengths  of 
between  10  and  50  km;  they  are  characterized  by  clockwise  rotations  when 
viewed  anti  parallel  to  F.  Curls  have  wavelengths  between  1  and  9  km  and  have 
lifetimes  between  0.4  and  2  sec.  Their  sense  of  rotation  Is  counterclockwise 
when  viewed  antiparallel  to  F.  Diameters  of  ~  2  km  are  typical  of  this  class 
of  vortices  (Webster  and  Hal 11nan»1973).  It  has  been  proposed  that  auroral 
curls  are  K-H  Instabilities  resulting  from  a  localized  excess  of  electrons 
within  a  sheet  of  precipitating  particles  (Hal  11  nan  and  Davis,  1970;  Webster 
and  Halil  nan,  1973).  The  two  vortex  structures  encountered  by  S3-2  most 
closely  resemble  auroral  curls. 

Webster  and  Hal  11  nan  (1973)  pointed  out  that  the  equations  for  K-H  insta¬ 
bilities  driven  by  either  unstable  current  sheets  or  unstable  charge  sheets 
have  the  same  form.  Three  aspects  of  our  observations  indicate  that  the 
vortices  reported  here  do  not  result  from  current  sheet  instabilities. 

First,  as  K-H  vortices  grow,  they  distort  the  east-west  alignment  of  the 
current  sheets.  The  S3-2  vector  magnetic  field  measurements  show  that  the 
intense  current  sheets  were  aligned  in  the  east-west  direction.  Second,  the 
vortices  are  detached  from  rather  than  embedded  In  the  current  sheets. 

Third,  the  sense  of  vorticlty  in  regions  where  electrons  carry  unstable 
currents  should  be  clockwise  (cf  Table  1  of  Webster  and  Hall i nan,  1973) 
rather  than  counter-clockwise  when  viewed  antiparallel  to  F  as  was  the  case 
In  our  observations. 

To  estimate  the  vorticlty  of  auroral  curls  HalHnan  and  Davis  (1970)  use 
the  growth  rate  for  K-H  Instabilities 
(3)  -  1  ^  [(1  -  2ka)2  -  e-4"*]1" 

where  k  Is  the  spatial  wave  number  of  the  curl  system  and  a  is  the  half-width  of 
the  auroral  charge  sheet.  Using  observatlonally  derived  values  u  ^  *  4.2  sec-1 
a  *  160  m  and  ka  *  0.15  in  equation  (3)  they  calculate  a  typical  vorticlty  of  28 
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4. 


sec"*.  This  Is  In  agreement  with  rocket  observations  of  vortlclty  (  u  ~  20  sec"*) 
near  the  edge  of  an  arc  (Kelley  and  Carlson,  1977).  Recently,  Wagner  et  al. 

(1981)  have  performed  computer  simulations  of  curl  formations  that  Include  the 
effects  of  partially  shielded  negative  charge  sheets.  In  this  case  the  K-H 
growth  rate  is 

(4)  »,-l  «,C(1  -  2  4ka)2  -  e-''*]1'2 

Equation  (4)  reduces  to  the  unshielded  case,  represented  by  equation  (3)  when 
the  parameters  a  *  1.  Their  simulations  show  that  for  increased  shielding 
(decreasing  A  )  the  maximum  growth  occurs  at  increased  values  of  ka  but  the  ratio 
u>  /  a)  s  decreases. 

Several  points  of  comparison  should  be  made.  First,  vortices  typical  of 
auroral  curls  calculated  by  Hallinan  and  Davis  greatly  exceed  those  estimated 
from  S3-2  measurements  of  the  auroral  electric  field  but  not  those  found  in 
rocket  measurements.  It  should  be  noted,  however,  that  the  auroral  sheet  thick¬ 
nesses  of  160  m  used  by  Hallinan  and  Davis  in  their  calculations  of  u  s  could 
not  be  resolved  at  sampling  rates  used  by  S3-2.  To  resolve  a  plasma  vortex  at 
least  two  electric  fields  samples  are  required.  For  S3-2  this  requires  1/16 
sec,  or  a  vortex  diameter  of  500  m.  The  vortices  detected  by  S3-2  had  radii  of 
2  km  and  0.75  km.  If  Hallinan  and  Davis  had  used  a  vortex  radius  of  1  km,  an 
averaqe  curl  wavelength  of  5  km  and  a  growth  rate  of  4.2  sec"*,  then  the  vorticity 
determined  from  equation  (3)  is  5.6  sec"*.  This  is  much  closer  to  our  estimated 
vortlcltles  of  1.2  and  2.5  sec"*.  Second,  In  the  simulations  of  Wagner  et 
al .  (1981)  and  in  the  rocket  measurements  the  curls  are  embedded  in  the  unstable 
charge  sheet.  We  have  already  pointed  out  that  the  vorticies  are  detached  from 
the  current  sheets  which  may  account  for  the  lower  vorticity.  Because  the  particle 
detector  was  sampling  low  energy  channels  as  S3-2  passed  through  the  vortices. 

It  Is  not  clear  whether  a  ground-based  television  camera  would  have  detected 
the  vortices  Inside  a  visible  auroral  arc. 


Another  point  of  Interest  Is  that  both  the  rocket  and  the  satellite  data 
Indicate  large-scale  electric  fields  In  the  range  of  100  mV/m.  This  is  con¬ 
siderably  less  than  the  estimate  of  500  mV/m  In  the  vicinity  of  auroral  curls  made 
by  Hallitnan  and  Davis  (1970).  The  most  likely  explanation  is  that  the  apparent 
motion  of  the  curls  Is  due  to  the  large  perpendicular  electric  fields  near 
the  source  of  auroral  beams  at  -  lRe  (Mozer  et  a!..  1977)  which  attenuate 
In  mapping  to  ionospheric  altitudes. 

Finally,  It  is  worth  reiterating  the  point  made  by  Kelley  and  Carlson  (1977) 
that  vortices  at  the  edge  of  auroral  arcs  may  create  some  portion  of  the  observed 
high-latitude  irregularities  in  plasma  density.  Since  it  is  likely  that 
horizontal  gradients  in  plasma  density  exist  near  auroral  arcs,  the  perturbed 
velocity  field  will  create  density  fluctuations  via  the  term 
an/  at  =  v  •  "v  n 

If  vortices  with  V  -  500  m/s  operate  on  a  10  km  density  gradient,  the  irregularity 
growth  rate  is  ,05  sec"*  which  is  higher  than  many  proposed  plasma  instabilities. 

In  the  introduction  of  this  paper,  we  noted  that  several  investigators 
have  studied  relationships  between  jt  and  *j.  A  field-aligned  potential 
drop  opens  the  atmospheric  loss  cone  for  individual  particles  in  ways  that 
depend  on  the  location  in  altitude  of  Assuming  an  isotropic,  Maxwellian 
parent  population,  Lyons  (1981)  has  shown  that  the  field-aligned  current 
density  is 

<5>  i  ,  '  J„  fi  [1-0-  |*>  e*P  1  -  14- }  3 

Bv  B1  kT(|i  -1) 

“V 

where  j0  is  defined  in  equation  (2);  B-j  and  Bv  are  the  magnetic  field  strengths 
at  the  lonsphere  and  at  the  top  of  the  voltage  drop,  respectively.  Information 
about  the  altitude  of  the  potential  drop  is  contained  in  the  Bj/Bv  ratio.  In 
the  limit  B-j /Bv  *  1,  the  potential  drop  is  located  just  above  the  ionosphere, 
and  j  |  «  J0  no  matter  how  large  the  potential  drop.  In  the  limit  q  /kT  •  , 


j  I*  j0  B^Bo.  Figure  10  Is  a  plot  of  j  (  /j0  versus  q  *|  /kl  for  six  values 
Of  B^/By. 

Crossing  the  morning  event  the  S3-2  magnetometer  measured  an  upward  current 
density  of  45  y  A/m2.  The  values  of  jQ  and  q  *,  /kT  Inferred  from  the  particle 
detector  measurements  are  1.2  y  A/m?  and  15,  respectively.  This  event 
designated  by  the  letter  M  In  Figure  10,  lies  well  above  all  of  the  current- 
voltage  characteristic  curves. 

In  the  evening  sector  event  we  measured  two  populations,  one  hot  (kT  = 

4  keV)  and  one  cold  (kT  *  200  eV).  We  now  attempt  to  evaluate  their  relative 
contributions  to  j  |  .  Magnetometer  data  presented  in  Figure  4  shows  that 
j  ,*  135  y  A/m2.  The  values  of  jQ  for  the  hot  and  cold  electrons  are  18.6  and 
3.45  yA/m2,  respectively;  ♦,  *  3.5  kV.  Using  the  symbols  H  and  C  in  Figure  9 
we  have  marked  the  current-voltage  relationship  one  gets  by  assuming  that 
either  the  hot  (H)  or  the  cold  (C)  electrons  Is  solely  responsible  for  j (  . 

Both  points  are  well  above  the  theoretical  curves.  Thus,  if  the  theoretical 
relationships  between  j  (  and  $|  are  valid  then  the  current  cannot  be  carried 
by  just  one  of  the  populations. 

We  next  assume  that  the  hot  population  carries  the  current  prescribed 
by  the  j  |  /  *!  relationship.  For  q  /kT  «  0.875,  j  ,  *  1.7  jQ  *  32  y  A/m2. 
Thus,  the  cold  population  must  be  responsible  for  103  y  A/m^.  This  is  plotted 
as  point  C'  in  Figure  10.  Point  C'  also  lies  well  above  any  characteristic  curve. 
Similarly  If  we  assume  that  the  cold  population  carries  the  current  prescribed 
by  the  j  i  /  *|  relationship,  we  find  that  the  hot  population  (point  H'  In 
Figure  10)  must  carry  almost  twice  the  predicted  current. 

In  order  to  reconcile  S3-2  measurements  with  theoretical  predictions  It  is 
necessary  to  assume  that  we  have  underestimated  either  or  n  In  analyzing  the 
electron  measurements.  For  example,  if  has  a  latitudinal  gradient,  a  conroon 
feature  of  Inverted-V  precipitation,  then  our  estimate  of  «,  would  only  be  a 
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lower  bound.  This  provides  a  viable  explanation  of  the  evening  side  event. 

Recall  that  the  evening  side  electron  burst  was  detected  0.25  sec  prior  to  (2  km 
poleward  of)  the  Intense  current  sheet.  He  find  that  If  rose  to  5.1  kV  above 
the  arc  then  the  observations  and  theory  agree.  The  points  plotted  as  H"  and  C“ 
show  that  with  *,  *  5. Ik V  the  hot  population  carries  a  current  of  45  u  A/ra‘  and 
the  cold  population  carries  100  u  A/m2. 

The  morning  side  electron  burst  exactly  coincides  with  the  45  p  A/m^ 
upward  current  sheet.  Thus,  a  gradient  in  $|  is  not  a  likely  explanation  of 
discrepancies  between  observations  and  theory.  However  an  underestimate  of 
the  parent  density  is  a  possible  explanation.  Recall  that  the  morning-side, 
burst  electrons  were  detected  while  the  detector  was  sampling  pitch  angles 
near  90°.  Burch  et  al.  (1979)  have  found  that,  at  least  in  the  case  of  polar 
cap  acceleration  events,  the  highest  levels  for  electrons  in  the  peak  channel 
are  found  along  the  direction  of  the  magnetic  field.  Thus,  using  a  value  f0 
measured  near  90°  pitch  angles  (Eq.(l))  to  calculate  n  could  lead  to  a  serious 
underestimate.  Agreement  with  theory  requires  that  the  parent  density  be 
0.75  cm“3  rather  than  our  estimated  value  of  0.3  cm~3.  In  both  the  morning 
and  evening  events,  the  theory  and  measurements  agree  only  if  the  ratio  B-j /8V 
is  relatively  high.  This  suggests  that  for  high  values  of  the  field- 
aligned  potential  drop  must  begin  at  high  altitudes. 

One  final  comment  is  appropriate.  Lyons  (1981)  stresses  the  fact  that 
collisionless,  single  particle  theory  takes  us  along  way  toward  understanding 
the  j  |  /  $|  relationship.  Conversely,  theories  of  anomalous  resistivity,  such 
as  that  proposed  by  Kindel  and  Kennel  (1971)  have  little  applicability  for 
auroral  currents.  In  the  case  of  upward  curents,  this  may  be  both  true  and  beside 
the  point.  The  model  of  Kindel  and  Kennel  concerns  plasmas  whose  ions  and 
electrons  have  a  relative  drift  motion  along  B,  not  a  beam  of  energetic  electrons 
passing  through  a  stationary  ionospheric  plasma.  The  model  does  apply  to  auroral 
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return  currents  which  are  mostly  carried  by  upwelllng,  cold  electrons.  Burke  et 
al .  (1980)  studied  the  case  of  a  10  w  A/m 2  return  current  from  a  discrete  arc. 

In  this  case  the  measured  plasma  parameters  satisfied  Kindel  and  Kennell  criteria 
for  marginal  stability.  Within  this  downward  current  sheet,  S3-2  detected  an 
electric  field  component  also  directed  along  B  into  the  ionosphere.  It  is  of 


interest  that  in  the  cases  presented  here  currents  out  of  the  ionosphere  were  45 
and  135  u  A/m^  while  those  directed  into  the  ionosphere  were  restricted  to  10  to 
15  v A/m^  range.  The  latter  numbers  require  cold  plasma  drifts  near  Kindel  and 
Kennel  limit  of  stability. 

Summary  and  Conclusions 

In  this  paper  we  have  presented  a  detailed  analysis  of  two  unusual  auroral 
events  observed  by  S3-2.  Both  events  were  characterized  by:  (1)  upward  current 
sheets  with  latitudinal  thicknesses  of  <  2  km  and  field-aligned  current  densities 
greatly  in  excess  of  a  few  y  A/m^,  (2)  return  currents  with  maximum  Intensities 

in  the  10  to  15  u A/m2  range,  with  the  principal  return  spatially  separated  for 
the  upward  current,  (3)  electric  field  rotations  driving  plasma  vortices  at 
locations  slightly  detached  from  the  upward  current  sheets,  and  (4)  bursts  of 
electrons  accelerated  through  potential  drops  of  several  kilovolts. 

The  counterclockwise  vortices  observed  by  S3-2  have  characteristics  similar 
to  auroral  curls.  The  results  are  consistent  with  theoretical  and  computer  models 
that  indicate  that  the  vortices  result  form  unstable,  auroral  charge  sheets. 

At  first  view  the  j  i  /  measurements  do  not  appear  to  agree  with  the 
predictions  of  collisionless,  single-particle  theory.  We  note,  however,  that 
If:  (1)  *|  had  a  maximum  value  of  5.1  kV  above  the  evening-side,  upward  current 
sheet,  and  (2)  the  density  of  the  morning  side  event's  parent  population  has 
been  underestimated,  then  agreement  between  theory  and  observations  is  possible. 

In  both  Instances  the  field-aligned  potential  must  extend  to  relatively  great 
altitudes  above  the  auroral  ionosphere. 


ACKNOWLEDGEMENTS 


We  wish  to  thank  R.  Bostrom  of  the  Uppsala  Ionospheric  Observatory, 

C.-G.  Falthammer  of  Royal  Institute  of  Technology,  Stockholm  and  M.C.  Kelley  of 
Cornell  University  for  helpful  discussions  and  comments  on  the  text.  This  work 
was  supported  in  part  under  AFOSR  contract  78-3731  with  Northeastern  University. 


References 


Alfven,  H.,  Cosmlcal  Electrodynamics,  Oxford  University  Press,  Oxford,  p.  206,  1950 

Burch,  J.L. ,  S.A.  Fields,  W.B.  Hanson,  R.A.  Heells,  R.A.  Hoffman  and  R.W.  Janetzke, 
Characteristics  of  auroral  electron  acceleration  regions  observed  by  Atmosphere 
Explorer  C,  J.  Geophys.  Res.,  81,  2223,  1976, 

Burch,  J.L.,  S.A.  Fields,  R.A.  Heells,  Polar  cap  electron  acceleration  regions, 

J.  Geophys.  Res.,  84  ,  5863,  1979. 

Burke,  W.J.,  Electric  Fields,  Birkeland  currents  and  electron  precipitation 
in  the  vicinity  of  discrete  auroral  arcs,  in  Physics  of  Auroral  Arc  Formation, 
ed.  by  S.-I.  Akasofu,  and  J.R.  Kan,  American  Geophysical  Union,  Geophysical 
Monograph  25,  Washington,  0.C.,  p.  164,  1981. 

Burke,  W.J.,  D.A.  Hardy,  F.J.  Rich,  M.C.  Kelley,  M.  Smiddy,  B.  Shuman,  R.C.  Sagalyn 
R.P.  Vancour,  P.J.L.  Wildman  and  S.T.  Lai,  Electrodynamci  structure  of  the 
late  evening  sector  of  the  auroral  zone,  J.  Geophys.  Res.,  85,  1179,  1980. 

Chen,  C.-K.,  R.A.  Wolf,  M.  Harel ,  J.L.  Karty,  Theoretical  magnetograms  based 
on  quantitative  simulation  of  a  magnetospheric  substorm,  submitted  to 
J.  Geophys.  Res.,  1982. 

Choy,  L.W.,  R.L.  Arnoldy,  W.  Potter,  P.  Kintner  and  L.J.  Cahill,  Jr.,  Field- 
aligned  currents  near  an  auroral  arc,  J.  Geophys.  Res.,  76,  8279,  1971. 

Davis,  T.N. ,  and  T.J.  Hallinan,  Auroral  spirals  1.  observations,  J.  Geophys. 

Res.,  81,  3  953,  1976. 

Frlden,  M. ,  and  J.  Lemalre,  Relationship  between  auroral  electron  fluxes  and 
field-aligned  electric  potential  difference,  J.  Geophys.  Res.,  85,  664,  1980. 

Goertz,  C.  K.  and  R.W.  Boswell,  Magnetosphere-ionosphere  coupling,  J.  Geophys. 

Res.,  84:  7239,  1979. 

Hallinan,  T.J.,  Auroral  spirals  2.  theory,  J.  Geophys.  Res.,  81 ,  3959,  1976. 

Hallinan,  T.J.,  and  T.W.  Davis,  Small-scale  auroral  arc  distortions.  Planet 


Harel,  M.,  R.A.  Wolf,  P.H.  Relff,  R.W.  Spiro,  W.J.  Burke,  F.J.  Rich,  and  M.  Smiddy, 
Quantitative  simulation  of  a  magnetospherlc  substorm,  1.  model  logic  and  overview 


J.  Geophys.  Res.,  86,  2217,  1981a. 

Harel,  M. ,  R.A.  Wolf,  R.W.  Spiro,  P.H.  Relff,  C-K.  Chen,  W.J.  Burke,  F.J.  Rich 
and  M.  Smiddy,  Quantitative  simulation  of  a  magnetospherlc  substorm, 

2.  comparison  with  observations,  J.  Geophys.  Res.,  86,  2242,  1981b. 

Karty,  J.L. ,  C.-K.  Chen,  R.A.  Wolf,  M.  Harel,  and  R.W.  Spiro,  Modeling  of  hiqh- 
latitude  currents  in  a  substorm,  J.  Geophys.  Res.,  87 ,  777,  1982. 

Kelley,  M.C.,  and  C.W.  Carlson,  Observations  of  intense  velocity  shear  and 
associated  electrostatic  waves  near  an  auroral  arc,  J.  Geophys.  Res.,  82 , 

2343,  1977. 

Klndel,  J.M. ,  and  C.F.  Kennell,  Topside  current  instabilities,  J.  Geophys.  Res., 

76,  3065,  1971 

Knight,  L. ,  Parallel  electric  fields.  Planet.  Space  Sci.,  21 ,  741,  1973. 

Lemaire,  J.  and  M.  Scherer,  Plasma  sheet  particle  precipitation:  a  kinetic  model. 
Planet.  Space  Sci , ,  21 ,  281,  1973. 

Lemaire,  J.  and  M.  Scherer,  Ionosphere-plasmasheet  field-aligned  currents  and 
parallel  electric  fields.  Planet.  Space  Sci.,  22,  1485,  1974. 

Lyons,  L.R.,  The  field-aligned  current  versus  electric  potential  relation  and 
auroral  electrodynamics,  in  Physics  of  Auroral  Arc  Formation,  ed.  by  S.-I.  Akasofu 
and  J.R.  Kan,  American  Geophysical  Union,  Geophysical  Monograph  25,  Washington 
D.C.,  p.  252,  1981. 

Lysak,  R.L.  and  C.W.  Carlson,  The  effect  of  microscopic  turbulence  on  magnetosphere 
ionosphere  coupling,  Geophys.  Res.  Lett.,  _8:  269,  1981. 

Lysak,  R.L,  and  C.T.  Dum,  Dynamics  of  magnetosphere-ionosphere  coupling  including 
turbulent  transport,  J.  Geophys.  Res.,  (accepted  for  publication,  1982). 


Mozer,  F.S.,  C.W.  Carlson,  M.K.  Hudson,  R.B.  Torbert,  B.  Parady,  J.  Yatteau, 
and  M.C.  Kelley,  Observations  of  paired  electrostatic  shocks  In  the  polar 
magnetosphere,  Phys.  Res.  Lett,,  38,  292,  1977. 

Mozer,  F.S.,  C.A.  Cattell,  M.K.  Hudson,  R.L.  Lysak,  M.  Temerln  and  R.B.  Torbert, 
Satellite  measurement  and  theories  of  low  altitude  auroral  particle  accele¬ 
ration,  Sp,  Scl.  Rev.,  27,  155,  1980. 

Ogutl,  T.,  TV  observations  of  auroral  arcs.  In  Physics  of  Auroral  Arc  Formation, 
ed  by  S.-I.  Akasofu  and  J.R.  Kan,  American  Geophysics  Union,  Monnograph  25, 
Washington,  D.C.  p.31,  1981. 

Sniddy,  M. ,  W.J.  Burke,  M.C.  Kelley,  H.A.  Saflekos,  M.S.  Gussenhoven,  D.A.  Hardy 
and  F.J.  Rich,  Effects  of  high  latitude  conductivity  on  observed  convection 
electric  fields  and  Birkeland  currents,  J.  Geophys.  Res.,  85,  6811,  1980. 

Spiro,  R.W.,  M.  Harel ,  R.A.  Wolf,  and  P.H.  Reiff,  Quantitative  simulation  of 
magnetospheric  substorm,  3.  plasmaspheric  electric  fields  and  evolution  of 
plasmapause,  J.  Geophys.  Res. ,  86,  2261,  1981. 

Wagner,  J.S.,  J.R.  Kan,  S.-I.  Akasofu,  T.  Tajima,  J.N.  Leboeuf  and  J.M.  Dawson, 

A  simulation  study  of  '/-potential  double  layers  and  auroral  arc  deformations, 
in  Physics  of  Auroral  Arc  Formation,  ed.  by  S.-I.  Akasofu  and  J.R.  Kan,  American 
Geophysical  Union,  Geophysical  Monograph  25  Washington,  D.C.,  p.  304,  1981. 

Webster,  H.F.  Structure  in  magnetically  confinded  electron  beams,  J.  Appl .  Phys. , 
28,  13  95,  1957. 

Webster,  H.F.,  and  T.J.  Hallinan,  Instabilities  in  charge  sheets  and  current 

sheets  and  their  possible  occurrence  in  the  aurora.  Radio  Science,  j8,  475,  1973. 

Whipple,  E.C. ,  The  signature  of  parallel  electric  fields  in  a  collisionless 
plasma,  J.  Geophys.  Res.,  82,  1525,  1977. 

Winningham,  J.O.,  F.  Yasuhara,  S.-I.  Akasofu  and  W.J.  Heikkila,  The  latitudinal 
morphology  of  10-eV  to  10-keV  electron  Fluxes  during  magnetically  gulet  and 
disturbed  times  in  the  2100-0300  MIT  sector,  J.  Geophys.  Res.,  80,  3148,  1975. 


Figure  Captions 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Sketch  of  satellite  pass  to  west  of  a  negatively  charged  flux  tube. 

The  electric  fields  and  T  x  "5  drifts  due  to  the  space  charge  electric 
field  Tc  and  the  main,  convection  electric  field  ?0  are  indidicated, 
the  panels  give  the  x  component  of  Ec  measured  by  the  satellite. 

The  potentials  due  to  the  space  charge  *c,  convection  $0  and 
their  superposition  4>s  as  functions  of  time  (distance)  along 
the  satellite  trajectory. 

Measurements  of  A  By  for  S3-2  passes  over  the  dawn-side,  southern 
hemisphere,  auroral  oval  during  Rev.  4079A/S,  (top)  and  Rev  4079  8/S 
(bottom).  Rough  estimates  of  region  1,  region  2  and  polar  cap 
positions  are  provided  for  reference. 

The  up-down  (UD),  north-south  (f6)  and  east -west  components  of  TT 
in  geomagnetic  coordinates  across  dawn-side  event.  Values  of  j  ( 
are  positive  for  currents  out  of  the  ionosphere.  Vector  representa¬ 
tions  of  the  electric  field  components  are  also  in  geomagnetic 
coordinates. 

Expanded  plot  of  morning  side  event.  The  top  three  plots  give  the 
east -west  magnetic  deflection,  the  north-south  electric  field  component 
and  the  potential  along  the  satellite  trajectory.  The  next  panel 
gives  the  T  x  f  drift  velocity  V.  The  bottom  panel  combines  the  ♦  and 
V  plots  to  approximate  the  two  dimensional  plasma  flow  in  the  vicinity 
of  the  upward  current  sheet. 

Five  consecutive  directional  differential  flux  spectra  observed  near  the 
morning-side  event.  The  UT's  at  which  spectra  accumulations  began  as 
well  as  the  pitch  angles  sampled  half  way  throuqh  the  one  second 
accumulation  periods  are  provided.  Oashed  lines  give  the  one  count 
per  accumulation-period  sensitivity  of  the  detector. 
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Figure  6.  Magnetic  field  deflections,  current  Intensities  and  electric  field 
vectors  for  evening-side  structure  In  the  same  format  as  Figure  3. 


Figure  7 


Figure  8 


Figure  9, 


Expanded  plot  of  the  potential  and  flow  pattern  In  the  vicinity  of 
the  evening-side  upward  current  sheet. 

Six  consecutive  spectra  measured  near  evening  side  event  in  the 
same  format  as  Figure  5. 

Semilog  plot  of  distribution  function  versus  energy  for  spectrum  3 
of  Figure  5  (top  panel)  and  spectrum  2  of  Figure  8  (bottom  panel). 
Dashed  lines  indicate  the  one  count  per  accumulation  period 
sensitivity  of  the  dectector. 


Figure  10.  Plot  of  j|  /j0  versus  kT/  for  various  Bv/Bi  ratios. 
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Abstract 

The  role  of  the  Pierce  instability  in  the  formation  of  double 
layers  is  considered  and  compared  with  that  of  the  Suneman  in¬ 
stability.  Pierce  instabilities  have  been  identified  in  a  double- 
-layer  experiment,  where  they  lead  to  ion  trapping.  Here  the 
effects  of  external  circuit  elements  are  considered.  In  the  case 
of  immobile  ions  the  onset  criteria  are  unaffected,  but  in  the 
unstable  range  the  growth  rate  is  reduced  by  the  external  impe¬ 
dance.  Required  experimental  values  of  the  circuit  elements  are 
estimated.  The  possible  relevance  to  computer  simulations  is 
noted. 
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1 .  Introduction 

Pierce  (1944)  showed  that  a  finite  length  system  consisting  of 
an  electron  beam  moving  through  a  neutralizing  background  of 
fixed  ions  can  be  made  unstable  if  the  end  grids  are  externally 
short  circuited.  If  the  ions  are  allowed  to  respond  the  insta¬ 
bility  still  occurs  (Faulkner  and  Ware,  1969)  ana  in  addition 
there  is  an  unstable  mode  depending  on  ion  oscillations  which 
is  dominant  for  long  systems  and  corresponds  to  the  ordinary 
two-stream  Buneman  instability.  The  effects  of  finite  ion  mass 
and  of  electron  thermal  velocity  spread  were  treated  by  Yuan 
(1977).  Saeki  et  al.  (1977)  analysed  the  transition  from  the 
Buneman  to  the  Pierce  mode  with  increasing  beam  current  and 
confirmed  their  results  experimentally.  From  a  nonlinear  treat¬ 
ment  of  the  Pierce  mode  Shapiro  and  Shevchenko  (1967)  find  that 
even  close  to  the  threshold  current  the  energy  of  the  disturbance 
becomes  comparable  with  the  initial  beam  energy. 

In  the  long-wavelength  regime  the  ordinary  two-stream  instability 
for  an  unbounded  plasma  can  drive  local  evacuations ,  which  may 
initiate  the  formation  of  double  layers  (Raadu  and  Carlqvist, 
1981).  The  local  plasma  density  in  the  evacuation  region  con¬ 
tinues  to  drop  even  in  the  non-linear  regime  as  shown  by  numeri¬ 
cal  analysis  including  the  ion  Vlasov  equation  (ibid)  and  in 
an  analytic  fluid  treatment  (Galeev  et  al. f  1981).  Quasineutrali¬ 
ty  holds  during  the  initial  evacuation,  so  that  the  ion  and 
electron  density  perturbations  are  approximately  equal.  The 
final  formation  of  a  double  layer  must  depend  on  departures 
from  charge  neutrality  and  possibly  on  particle  trapping  as 
indicated  by  numerical  solutions  of  the  complete  Vlasov-Poisson 
system  of  equations  (Chanteur  and  Volokhitin,  1982). 

Iizuka  et  al.  (1979)  have  observed  the  formation  of  a  double 
layer  as  the  result  of  the  injection  of  an  electron  beam  into 
a  low  density  plasma.  As  they  also  described  in  the  preliminary 
report  (Saeki  et  al. ,  1977)  Pierce  instabilities  supplant  the 
Buneman  instability  beyond  a  critical  current  level.  They 
ascribe  the  formation  of  a  potential  minimum  to  the  rapid  growth 
of  the  Pierce  mode.  Ions  are  trapped  in  this  suddenly  formed 
potential  well  and  constitute  the  thermal  population  on  the 
low-potential  side  of  the  double  layer.  Thus  in  their  experi¬ 
ment  the  Pierce  instability  plays  a  vital  role  in  the  formation 
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of  the  double  layer. 

The  essential  difference  between  the  situation  considered 
by  Raadu  and  Carlqvist  (1981)  and  that  investigated  by 
Iizuka  et  al .  (1979)  is  the  length  of  the  system.  For  a  long 
system  mam y  unstable  Buneraan  inodes  are  present  and  the  form  of 
<&•  growing  perturbations  reflects  their  initial  random  nature, 
density  depletions  leading  to  local  evacuations  which  may  ini¬ 
tiate  double-layer  formation.  In  a  short  system  a  single  Pierce 
mode  can  dominate  and  determine  the  form  of  the  growing  pertur¬ 
bation  and  in  particular  the  tendency  to  form  local  evacuation 
regions  is  lost. 

In  view  of  these  considerations  it  is  of  interest  from  the  point 
of  view  of  double-layer  formation  to  look  further  into  the 
conditions  for  the  growth  of  the  Pierce  instability.  Here  we  will 
investigate  the  influence  of  an  external  circuit  with  finite 
impedance  instead  of  the  normally  assumed  zero  impedance 
short-circuit.  In  a  laboratory  experiment  it  is  possible  to 
introduce  external  circuit  elements  and  to  influence  the 
growth  of  Pierce  instabilities.  In  a  cosmic  plasma  an  external 
current  system  feeding  a  region  of  double-layer  formation  may 
be  represented  by  an  equivalent  circuit.  In  both  cases  the  Pierce 
instability  may  be  suppressed  or  at  least  the  growth  rate  may 
become  insignificant . 

2.  Influence  of  an  external  circuit 

In  order  to  highlight  the  effects  of  an  external  circuit  we  will 
here  only  consider  the  case  of  a  cold  electron  beam  and  immobile 
ions.  For  an  unbounded  plasma  there  are  only  stable  Langmuir 
modes  given  by  the  standard  dispersion  relation. 


(w-Jcv)  2 

where  k  and  o  are  the  wave  number  and  frequency,  v  is  the  beam 
velocity  and  a>  the  plasma  frequency.  For  a  finite  length  sys¬ 
tem  with  space  coordinate  x  potential  perturbations  of  the  form 


$(x)  ■  [A  exp(ik+x)  +  E  exp(i)^x)  -  EQx)  e 


-iiat 
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become  possible  where  the  uniform  electric  field  EQis  due  to 
charges  on  the  end  grids  (x  =  o,  d)  and  k,,.  are  the  solutions 
of  the  free-mode  equation  (1). 

At  the  boundary  where  the  electron  beam  is  injected  (x  *  0)  the 
density  and  velocity  perturbation  are  set  to  zero.  Alternative 
boundary  conditions  {e.g.  setting  the  electrical  field  to  zero) 
are  in  general  not  acceptable,  since  they  imply  a  non-physical 
coupling  between  the  conditions  within  the  plasma  region  and 
the  parameters  of  the  beam  injection  mechanism.  The  electron 
beam  is  disturbed  by  the  external  circuit  only  through  the 
uniform  electric  field  EQ.  The  external  circuit  which  is  connected 
to  the  end  grids  responds  only  to  the  potential  difference 
♦  (d)  -  4(0)  .  There  is  no  direct  coupling  between  the  external 
current  which  builds  up  charges  on  the  end  grids  and  the  current 
carried  by  the  electron  beam.  This  severely  restricts  the 
possible  applications  to  cosmic  situations  where  in  general  the 
electron  current  within  a  region  of  double-layer  formation  is, 
together  with  the  external  current,  part  of  the  same  extended 
current  system  and  "end  grids"  which  may  carry  charge  do  not 
exist. 

In  the  standard  Pierce  analysis  the  two  boundary  conditions  with 
the  short-circuit  condition  4(d)  =  ♦(0)  lead  to  a  characteristic 
equation  determining  the  growth  rate  for  given  plasma  parameters 
(Pierce,  1944;  see  also  Mikhailovskii ,  1974).  Here  we  instead 
consider  the  external  circuit  illustrated  in  Figure  1 ,  consisting 
of  capacitance  C,  resistance  R  and  inductance  L  in  series.  For 
an  external  current  1=6  where  iQ  are  the  charges  on  the  end 
grids,  the  potential  across  the  plasma  satisfies 

4(d)  -4(0)  =  —  +  IR  +  L  —  (3) 

C  dt 

Introducing  the  vacuum  capacitance  of  the  plasma  volume 
Cc  **  cQA/d  where  A  is  the  area  of  a  grid  we  have  from  (2) 

A[exp(ik+d)-11  ♦  B[exp(ik_d) -1  ] 

«-2-  +  -2-+IR+L^  (4) 

C0  C  dt 


In  this  form  the  capacitance  Cc  enters  as  a  circuit  element  in 
series  with  those  ir.  the  external  circuit  and  is  thus  signifi  ¬ 
cant  for  determining  their  relative  importance .  Equation  (4) 
together  with  the  standard  bounder-  conditions  yields  the 
modified  characteristic  equation 


F(0f)  =  2  (.* 9  *  [1  -  e  6  cos  0g  ]  +  (-3^  -  fl  ')  e  9 sin  ae 

+  [1  +  c"1  t  SO  +  L32]  62  (3  2  +  0  2)  =  0  (5) 

e 


where  we  use  Yuan's  (1977)  notation  (3  - 
and  have  non-.iimensionabized  the  circuit 
R  *=  R  CQv/d ,  L  =  LCov2/d2).  If  C~1  *  R 
standard  characteristic  equation. 


-isod/'v,  8  =  id  d/v  ) 

parameters  (C  *  C/C0, 
=  L  *  0  we  recover  the 


3.  Discussion 

We  first  of  all  notice  that  foi*  marginal  instability  (3  =  0) 

it  follows  from  Equation  (5)  that  9e  «  mr.  This  condition  is 

identical  to  the  standard  case.  The  marginally  unstable  modes 

consist  of  an  integral  number  of  half  wave  lengths  of  spacially 

undamped  sinusoidal  oscillations  with  no  uniform  field  component 

(E  *  0) .  Clearly  the  external  circuit  cannot  influence  these 
^  » 
modes  as  there  is  no  required  external  current  (Q  =  0) . 

In  the  unstable  parameter  ranges  the  growth  rate  a  is  reduced 
by  the  external  impedances.  Thus  in  Figure  2  the  effect  of  an 
external  inductance  is  shown.  For  the  particular  choice  L  *  40 
the  maximum  growth  rate  is  reduced  by  a  factor  a0.53  as  compared 
with  the  case  of  an  external  short-circuit. 

The  growth  rate  9  decreases  rather  slowly  as  the  external  induc¬ 
tance  is  increased.  This  may  be  seen  from  Table  1,  where  the 
equivalent  necessary  values  of  the  inductance  L,  resistance  R 
or  capacitance  C  are  given  for  a  number  of  growth  rates 
(8  =  1.6u).  In  Equation  (5)  the  circuit  parameters  appear 

additively  in  the  3ame  term  and  hence  for  a  given  growth  rate 
it  is  easy  to  find  parameter  values  which  give  the  same  result. 
For  large  values  of  the  impedance-  we  have  the  approximations  , 


8  a  (a/L)~V*  =  (a/R)"V»  =  (aC)'1/1 


where  a  =  -  9  sine  ,  which  aciain  illustrates  the  relatively 
weak  dependence  on  the  external  circuit  parameters. 

For  the  experiment  of  Saeki  et  ad.  (1377;  lizuka  et  al . ,  1979) 
we  estimate  the  vacuum  capacitance  of  the  plasma  volume  to  be 
C0  a  0.25  pF  and  the  transit  time  d/v  =  10“^s.  To  determine 
the  size  of  external  circuit  elements,  which  according  to  the 
theory  presented  here  could  have  an  effect  if  introduced  into 
the  experiments,  we  note  that  for  the  dimensionless  parameters 
L,  R  and  C  to  be  equal  to  unity  the  actual  values  should  be 
L  -  0.04  H,  R  =  0.4J$Ohm,  C=  0.25  pF  .  In  particular  for  the 
theoretical  curve  giver,  in  Figure  2  (L  *  40)  the  corresponding 
outer  inductance  would  be  L  *  1.6  K.  The  required  dimensions  of 
the  circuit  elements  are  such  that  unless  they  are  deliberately 
introduced  into  the  experiment  the  impedance  of  the  outer  circui 
has  almost  certainly  a  negligible  effect  in  the  experiments  of 
Saeki  et  al.  (ibid) . 

In  view  of  the  comments  in  Section  2  the  special  configuration 
required  for  the  Pierce  instability  is  unlikely  to  be  applicable 
in  cosmic  situations.  However  if  we  estimate  the  external 
inductance  by  where  l  is  the  length  of  the  currant  system, 
and  consider  a  local  region  of  length  d  and  cross  section  d2 
the  dimensionless  inductance  L  =  (v/c)  2 (i/d).  Since  v  <  c  and 
we  should  have  i  >  d  there  is  no  a  priori  restriction  on  L. 

From  the  analysis  presented  here  we  must  have  h  >>  1  for  a 
significant  reduction  of  the  Pierce  instability  growth  rate. 

Finally  we  note  that  computer  simulations  of  double-layer  forma¬ 
tion  are  of  necessity  made  for  a  finite  length  system.  Instabili 
ties  of  the  Pierce  type  are  then  possible  and  may  than  be 
suppressed  by  including  the  effects  of  a  suitable  external 
circuit  in  the  simulation.  A  modified  form  of  the  Pierce  insta¬ 
bility  might  be  expected  for  example  if  a  constant  voltage  is 
imposed  across  the  simulation  region  (cf.  Smirnov,  1966).  This 
implies  that  the  perturbations  in  the  total  potential  difference 
must  be  zero  and  an  appropriate  uniform  electric  field  must  be 
suoerimoosed  as  in  the  standard  Pierce  instability  analysis. 
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Fig.  1 .  The  physical  system.  An  electron  beam,  (e  )  traverses  a 
region  or  stationary  ions  bounded  by  grids  (located  at 
x  »  0,  d)  carrying  charges  *Q.  The  charges  are  supplied 
via  an  external  circuit  with  current  I  passing  capaci¬ 
tance  C,  resistance  R  and  inductance  L  in  series. 

Fig.  2.  The  dimensionless  growth  rate  e(=  -io»d/v)  is  given  as  a 
function  of  the  electron-beam  parameter  (*  <i)ped/v} 
for  the  case  of  an  external  dimensionless  inductance  L  * 
*  40,  solid  curve,  and  for  comparison  for  the  case  of  an 
external  short-circuit  (L  =  0} ,  dashed  curve. 
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Abstract 


It  is  known  than  large  amplitude  oscillations  can  occur  in  the 
current  flowing  through  a  plasma  diode,  typically  when  a  constant 
potential  is  applied  across  the  device.  Burger  (1965)  suggested 
via  a  computer  simulation  that  the  oscillation  characteristics 
was  a  function  of  the  quantities  and  t.  ,  namely  the  respective 
time  for  an  electron  and  an  ion  to"  crosi  the  electric  field  re¬ 
gion  inside  the  diode.  On  the  rapid  time  scale  t  ,  the  self  con¬ 
sistent  equilibrium,  configuration,  was  unstable . eNorris  (  1  964) 
had  previously  arrived  at  the  same  conclusion  using  analytical 
arguments.  In  that  work,  it  was  concluded  that  the  instability 
occurred  since  the  diode  acted  as  a  negative  resistance  on  the  T 
scale.  A  positive  feedback  effect  forced  the  system  away  from  e 
equilibrium. 

Silevitch  (1981)  used  the  Burger  mechanism  to  suggest  an  expla¬ 
nation  for  the  flickering  aurora  phenomenon.  He  extended  the 
Norris  argument  and  showed  by  a  variational  method  that  a  plau¬ 
sible  analytic  model  for  a  double  layer  (DL)  behaved  as  a  nega¬ 
tive  resistance  on  the  scale.  In  this  present  work  we  re-exa- 
raine  the  negative  resistance  calculation  by  taking  a  more  detailed 
account  of  the  constraints  which  are  imposed  on  the  electron 
distributions  that  exist  in  the  DL  region.  Specifically ,  we  shall 
focus  upon  the  role  of  the  energetic  trapped  electrons  which  ori¬ 
ginate  at  the  high  potential  side  of  the  DL. 


♦Permanent  address:  Dept  of  Electrical  Engineering, 
Northeastern  University,  Boston,  MA  02115,  U.S.a! 


is  indeed  proportional  to  the  time  for  an  ion  to  traverse  a 
distance  which  is  characterized  by  the  spatial  extent  of  the 
entire  system  (i.e.  the  electrode  separation  distance).  This  is 
in  contrast  to  the  suggestion  of  Silevitch  (1981)  (henceforth 
referred  to  as  Paper  1)  that  the  transit  distance  comprise  only 
the  narrow  DL  region. 

In  this  presenc  work  we  shall  focus  upon  the  dynamics  of  the  DL 

on  the  r  time  scale.  We  will  reexamine  the  computation  of  ne- 
e  * 

gative  resistance  given  An  Paper  1  by  taking  a  more  detailed 
account  of  the  constraints  which  are  imposed  on  the  various 
electron  distributions  which  exist  within  the  DL  region.  It 
should  also  be  emphasized  that  the  DL  instability  discussed  here 
is  a  global  rather  than  a  local  effect.  By  this  we  mean  that 
the  properties  of  the  external  circuit  surrounding  the  DL  ele¬ 
ment  play  a  crucial  role  in  the  devlopment  of  the  relaxation 
oscillations.  Thus  in  a  given  application  it  is  not  enough  to 
determine  that  the  DL  behaves  as  a  negative  resistor  on  the  t 
scale.  One  must  also  examine  how  this  resistance  couples  through 
an  external  circuit  and  whether  it  causes  a  positive  feedback 
effect. 

As  in  Paper  1  we  shall  choose  a  particular  steady  state  DL  model 
and  demonstrate  that  it  exhibits  a  negative  resistance  on  the  re 
time  scale.  In  our  case  we  continue  using  the  Kan  and  Lee  (1980) 
model  introduced  in  Paper  1 .  It  is  intended  that  the  calculation 
presented  here  be  used  as  a  guide  if  one  wants  to  investigate 
the  stability  of  a  different  DL  model  that  may  better  describe 
some  particular  experimental  situation. 

A  similar  but  somewhat  abbreviated  version  of  this  paper,  Raadu 
and  Silevitch  (1982),  will  appear  as  part  of  the  proceedings  of 
a  small  symposium.  We  are  presenting  the  results  again  here 
because  the  earlier  version  will  not  be  widely  circulated  within 
the  plasma  physics  community.  Moreover  we  will  make  extensive 
use  of  the  material  already  contained  in  Paper  1 ,  in  order  to 
avoid  unneccessary  mathematical  detail. 

In  Section  2  we  present  a  summary  of  the  Kan  and  Lee  (1980)  DL 
model.  We  examine  the  perturbation  of  the  DL  on  the  te  time 
scale  in  Section  3.  The  variational  technique  given  in  Paper  1 
is  generalized  and  we  compute  the  negative  resistance  of  the  DL 
which  results  when  various  constraints  are  imposed  on  the 
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electron  distribution  functions.  Ic  is  also  found  that  the  immo¬ 
bile  ions  introduce  an  additional  term  in  the  perturbation  ex¬ 
pansion  that  was  overlooked  in  Paper  "I .  We  investigate  this 
term  in  Section  4  and  present  arguments  to  show  that  it  has 
a  negligible  effect  on  the  results  presented  here  and  in  Paper  1 . 


2 .  Equilibrium  model 

The  equilibrium  auroral  DL  model  used  in  Paper  1  is  taken  from 
the  work  of  Kan  and  Lee  (1980)  .  The  DL  potential 
structure  <Mx)  ,  varies  from  $  (x=0)  *  0  to  <p  (x=d)  =  $Q.  Moreover, 
as  shown  in  Figure  4  of  Paper  1 ,  there  are  three  distinct  electron 
populations  associated  with  the  Kan  and  Lee  DL  model.  These  are: 

(a)  Streamin2_electrons_qriginating_at_ j> j(x=0 


A  waterbag  velocity  distribution  is  chosen  for  this  population 
which  is  simply 
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Here  N 


where  V^(4)  =  (2 | e ■ 4/m) ^  and  Vu(4)  =  [4V0^  +  ( <t> )  1  ^  •  ucis  l,el 

and  V  j  are  the  density  and  streaming  velocity  of  the  electrons 
originating  at  x  =  0.  Following  Paper  1  we  can  easily  obtain 
expressions  for  the  electron  density,  n0^($)  and  its  integral 

at  any  point  4(x)  within  the  DL  region. 
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Note  that  it  would  be  more  precise  to  indicate  these  quantities 


as  functionals  of  f  . (N  ,  , 
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(b)  Tragged_degraded_grtmarY_and_secondary_electrons_originating 

at  4  (x*d)  =  4 
- - -o 

Again  for  these  particles  we  use  a  waterbag  distribution  f  . 
of  value  ffc  centered  at  v  =  0  and  cut  off  at  v  =  ±(2jej$/m)J.  Thus, 
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(2) 


where  N  is  the  density  of  these  trapped  electrons  at  x  -  d. 
Again  simple  expressions  for  nej-(4)  and  g0^  (4)  can  be  ob¬ 
tained  (see  Paper  1). 

(c)  Tragged_  low_energY_e  lectrons  or iginat  ing_at  _£j_x*d)__*  _ 

A  Maxwellian  distribution  is  assumed  for  this  population.  It  is 

a  full  range  function  characterized  by  the  parameters  N  f>>N  ,) 

eo  e  1 

and  kT0O (<< j e j 4Q)  which  respectively  represent  the  electron  den- 
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sity  and  thermal  energy  at  x  =  d. 


3 .  Perturbation  Analysis 


If  the  equilibrium  DL  structure  is  perturbed  on  the  te  time 
scale  then  the  dynamic  resistance  of  the  DL  is  defined  as 


R0 


(3) 


Here  A  is  the  DL  cross  sectional  area  and  <5j  (3$o)  represents 
the  electron. current  density  (potential)  variations  from  the 
equilibrium  values 

[i.e.  $  ( x )  **{x)  *  <Mx)  +  <5$(x)  ] .  Cleary  only  the 

electrons  in  category  (a)  contribute  to  j  and  thus, 

Sj  *  fiUe|NelVel)  *  ie|Vel6Nel  -  |e|Nel5Vel  (4) 

In  order  to  calculate  Rq  we  need  the  key  result 

0  -  «(9el  *  9et  +  9e0>  (5> 


This  equation  is  obtained  by  first  multiplying  Poisson's  equa¬ 
tion  by  d<j>/dx  and  then  integrating  from  x  =  0  to  x  =  d  assuming 
charge  neutrality  at  both  endpoints.  Finally,  the  same  procedure 
is  repeated  for  the  perturbed  state  (4*)  and  the  two  equations 
are  subtracted  keeping  only  first  order  terms  in  the  variations 
like  6$.  It  should  be  noted  that  to  this  order  in  5<t>  it  is  not 
neccessary  to  impose  strict  charge  neutrality  for  the  perturbed 
state  at  x  =  0,  d.  Moreover,  a  rigorous  derivation  of  (5)  would 

include  on  the  rhs  a  static  ion  term  <5G.  defined  by 

ion 


where  n^tx)  is  the  ion  density  profile  in  the  unperturbed  state. 
In  Paper  1  this  term  was  neglected.  An  argument  justifying  this 
approximation  is  presented  in  the  final  portion  of  this  paper. 

In  order  to  obtain  an  expression  for  Rp  we  need  to  specify  in 
detail  those  constraints  that  apply  during  the  initial  disrup¬ 
tion,  To  illustrate  this  consider  electron  population  (a) . 

Eq.  (4)  defines  a  reJation  between  <5Ne^  and  6Ve^.  Another  is 
needed.  For  example,  in  Paper  1  it  was  assumed  that  and 
V  .  were  independent  and  so  either  6N  .  *  0  (i.e.  strict  charge 
neutrality)  or  SV^,  =  0.  Perhaps  a  more  realistic  constraint 
would  be  to  impose  the  condition  5f1  *  0.  This  would  imply  from 
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Eq .  ( 1 ) 
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From  Eqs.  (4)  and  (7)  we  find 
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The  expression  for  R^  is  now  obtained  by  expanding  Eq.  (5)  as: 
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Note  that  these  constraints 

would 

ww 

not  impose  strict  charge  neutrality  at  x  =  d.  Using  Eqs.  (4)  and 
(9)  we  obtain  the  resistance  given  by 

(10) 
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If  we  replace  SN^  *  0  by  condition  (7)  we  find  the  DL  negative 

(11) 


resistance  R2  will  have  the  value  R2  where 


**  - ^  ft 


For  the  parameters  of  the  auroral  example  in  Paper  1  we  find 
that  the  two  terms  in  the  bracket  have  roughly  the  same  magni¬ 


tude  and  hence  R2  M  . 


Let  us  now  generalize  the  constraints  on  5N  .  and  5N  .  An  ob- 
3  et  eo 

vious  first  step 
This  would  imply 


vious  first  step  is  to  assume  strict  charge  neutrality  at  <fr  *  $Q. 
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For  simplicity  let  U3  also  still  assume  <5N  ,  =  5N 

el  eo 

Eqs.  (4)  and  (12)  we  find  that 
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the  DL  negative  resistance  will  now  have  the  value  R-j  given  by: 
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For  the  parameters  of  the  auroral  DL  of  Paper  1  we  find  that 
only  the  An^/SV^  correction  term  is  important  and 
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Moreover,  if  we  assumed  that  is  was  primarily  the  low  energy 

population  which  reacted  to  maintain  charge  neutrality  (i.e. 

6Nefc  =  0,  6Ngo  =-5ng^  (^q)  )  ,  then  the  resulting  negative  resistance 

would  be  approximately  equal  to  R1 .  A  result  similar  to  Eq.  (15) 

is  obtained  if  we  relax  the  charge  neutrality  constraint'  and 

assume  instead  Sf.  =  AN  -  0. 

t  eo 

From  the  above  discussion  we  conclude  that  the  response  of  the 
trapped  particle  populations  will  have  quite  an  important  effect 
upon  the  value  of  DL  negative  resistance.  According  to  the 
theory  in  Paper  1  a  smaller  value  of  negative  resistance  could 
quench  the  DL  disruption  on  the  Tg  time  scale.  To  test  this 
hypothesis  we  can  envisage  an  experiment  which  allows  the 
controlled  injection  of  trapped  electrons.  One  could  then  study 
the  DL  disruption  characteristics  as  a  function  of  the  trapped 
electron  distributions. 


4.  Approximate  Treatment  of  Unperturbed  Ions 

Before  concluding  this  paper  we  will  justify  the  neglect  of 
AG^Eq.  (6))  in  Eq.  (c).  The  unperturbed  electric  field  of  a 
strong  DL  structure  is  primarily  nonzero  in  an  interior  region 
of  space  Ax(<d) .  Under  the  assumption  that  the  perturbed  electric 
field  (— g~)  is  also  confined  to  essentially  the  same  Ax  region 
we  can  rewrite  Eq.  (6)  as, 

<6.  -  ^  fo  dx  U*  -  4)  *  n.  «4q. 

Here  n^  is  a  mean  value  of  unperturbed  ion  density  in  the  Ax  region 
It  is  given  by  the  expression. 


n1 


± 


l<*o> 


For  the  auroral  DL  model  discussed  in  Paper  1 ,  it  can  be  shown 
that  —  0.1  (3ge/34Q)  and  so  for  this  case  can  indeed 
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be  neglected  in  Eq.  (5). 
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Abstract 

It  is  known  that  large  amplitude  oscillations  can  occur  in  the 
current  flowing  through  a  plasma  diode,  typically  when  a  constant 
potential  is  applied  across  the  device.  Burger  (1965)  suggested 
via  a  computer  simulation  that  the  oscillation  characteristics 
was  a  function  of  the  quantities  re  and  r namely  the  respective 
time  for  an  electron  and  an  ion  to  cross  the  electric  field  re¬ 
gion  inside  the  diode.  On  the  rapid  time  scale  -tg,  the  self  con¬ 
sistent  equilibrium  configuration,  was  unstable.  Norris  (1964) 
had  previously  arrived  at  the  same  conclusion  using  analytical 
arguments.  In  that  work,  it  was  concluded  that  the  instability 
occurred  since  the  dicde  acted  as  a  negative  resistance  on  the 
scale.  A  positive  feedback  effect  forced  the  system  away  from 
equilibrium. 

Silevitch  (1981)  used  the  Burger  mechanism  to  suggest  an  expla¬ 
nation  for  the  flickering  aurora  phenomenon.  He  extended  the 
Norris  argument  and  showed  by  a  variational  method  that  a  plau¬ 
sible  analytic  model  for  a  double  layer  (DL)  behaved  as  a  nega¬ 
tive  resistance  on  the  scale.  In  this  present  work  we  re-exa¬ 
mine  the  negative  resistance  calculation  by  taking  a  more  detailed 
account  of  the  constraints  which  are  imposed  on  the  electron 
distributions  that  exist  in  the  DL  region.  Specifically,  we  shall 
focus  upon  the  role  of  the  energetic  trapped  electrons  which  ori¬ 
ginate  at  the  high  potential  side  of  the  DL. 


Key  words:  Double  Layers,  Stability,  Negative  Resistance, 
Relaxation  Oscillations,  Stationary  Ions. 
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Abstract 


Second  order  effects  are  calculated  for  a  low  frequency 
electromagnetic  instability  due  to  an  ion  beam  in  a  plasma. 
This  instability  is  a  characteristic  of  a  parallel  shock 
model  which  includes  an  ion  reflecting  electrostatic  sub¬ 
shock.  The  analysis  is  similar  to  that  done  by  S.  P.  Gary 
but  a  counterstreaming  configuration  is  chosen  which  is  homo¬ 
geneous  in  time  and  which  has  spacially  growing  modes.  Among 
other  effects,  energy  is  transf erred  from  the  incoming  main 
plasma  ions  to  the  beam  ions.  Electron  pressure  effects 
are  calculated  and  are  small  for  low  frequencies.  The 
application  of  this  model  to  the  earth's  bow  shock  is  dis¬ 
cussed. 
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1 .  Introduction 


A  shock  model  which  consists  of  an  ion  reflecting  potential 
barrier  immersed  in  electromagnetic  turbulence  has  been  sug¬ 
gested  in  theory  [Jackson  and  Golden,  1979]  and  by  computer 
simulation  of  oblique  shocks  [Lindman  and  Drummond,  1971; 

Biskamp  and  Welter,  1972].  The  simulation  results  of  Lindman 
and  Drummond  indicate  the  formation  of  an  ion  reflecting  po¬ 
tential  barrier  which  they  associate  with  the  ion  acoustic  wave. 
This  barrier  reflects  ions  upstream  causing  a  streaming  electro¬ 
magnetic  instability  which  is  identified  by  the  authors  as  a 
whistler  instability  (u>  >>  .  In  a  later  simulation,  Biskamp 

and  Welter  note  that  this  whistler  instability  stabilizes  by 
rearranging  the  electron  distribution  and  not  the  ion  distribu¬ 
tion.  They  suggest  that  the  ions  are  thermalized  by  the  occur¬ 
rence  of  a  nonlinear  two-stream  ion  acoustic  instability  which 
could  be  triggered  by  potential  variations.  This  instability 
would  then  heat  the  ions.  As  the  oblique  shock  becomes  more 
parallel,  the  high  frequency  whistler  instability  used  by  these 
authors  becomes  stable.  It  was  later  shown  [Jackson  and  Golden, 
1979]  that  for  parallel  shocks  the  low  frequency  whistler 
(u>  <  streaming  instability  takes  the  place  of  the  high  fre¬ 
quency  whistlers  used  in  the  oblique  case. 

We  note  further  that  the  reflecting  potential  model  extends 
the  high  beta  parallel  shock  model  of  Kennel  and  Sagdeev  [1967] 
to  the  lower  beta  (-0(1))  range.  In  the  Kennel -Sagdeev  model  a 
relatively  long  scale  length  potential  barrier  causes  an  ion 
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pressure  anisotropy  which  drives  Alfvdn  waves  unstable.  For 

lower  beta,  the  firehose  instability  is  no  longer  dominant.  In 

addition,  the  potential  barrier  becomes  more  narrow  and  greater 

in  magnitude  due  to  the  increase  in  shock  speed  relative  to  the 

2 

ion  acoustic  speed  (Cg  =  <  Te/m^) .  It  is  suggested  here  that 
this  larger  potential  reflects  ions  into  the  upstream  region 
and  that  the  firehose,  fluid,  Alfv6n  instability  is  replaced  by 
a  two-stream,  low  frequency,  whistler  instability. 

This  paper  is  an  examination  of  the  second  order  effects 
of  the  low  frequency  turbulence  in  the  reflected  ion  shock  model. 
It  is  similar  to  the  work  of  Barnes  [1970]  but  the  energy  and 
momentum  transfers  between  each  fluid  are  examined  in  more  detail 
and  the  second  order  electric  field  is  retained.  First  we  review 
the  notation  and  the  linear  theory.  Then  the  second  order  ef¬ 
fects  are  calculated  for  spacially  growing  modes  having  time 
independent  average  properties.  The  effect  of  this  turbulence 
on  the  pre-barrier  plasma  is  analyzed  using  the  second  order 
theory  of  Gary  and  Feldman  [1978;  Gary,  1978].  Their  theory  for 
temporal  growth  is  modified  here  for  the  case  of  spacially 
growing  modes  and  warm  electrons.  Also,  the  expressions  derived 
do  not  require  the  assumption  u  <<  which  was  used  by  Gary 
[1978].  Finally,  the  similarities  between  this  model  and  ob¬ 
servations  of  the  earth's  bow  shock  are  noted. 

2.  Linear  Analysis 

In  our  previous  work  we  supposed  the  existence  of  a  poten¬ 
tial  barrier  which  reflects  a  small  number  of  ions  into  the 
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pre-barrier  region.  These  reflected  ions  counterstream  with  the 
incoming  plasma  ions  causing  a  resonant  low  frequency  whistler 
instability.  The  linear  stability  analysis  assumed  a  locally 
homogeneous  plasma  where  modes  may  grow  in  time.  Low  frequency 
whistlers  become  unstable  for  Alfv^n  Mach  numbers  greater  than 
about  1.5  if  it  is  assumed  that:  (i)  the  ion  acoustic  soliton 
is  the  reflecting  potential,  (ii)  the  Alfv6n  and  ion  acoustic 
speeds  are  equal,  and  (iii)  electron  trapping  in  the  potential 
is  allowed. 

It  is  assumed  that  this  very  steep  potential  barrier 

exists  at  a  point  2=0  (see  Figure  1)  and  reflects  a  small 

amount  of  the  cool  incoming  plasma  into  the  upstream  region 

(2  <  0)  along  a  constant  magnetic  field,  Bq  =  Bq  2.  At  the 

point  z  -  0~  the  reflecting  potential  is  assumed  to  be  rero. 

The  ion  distribution  which  results  would  then  be  similar  to 

1/2 

Figure  2  where  v  =  (2e  ♦_/m4 )  '  and  *  is  the  maximum  poten- 
3  m  m  1  Tm 

tial  of  the  barrier. 

Using  this  ion  distribution  in  the  kinetic  dispersion 
relation  for  k| |Bq,  we  make  the  following  assumptions:  (i)  the 
incoming  ion  stream  and  the  electrons  appear  cold  to  the  un¬ 
stable  waves,  (ii)  |u|  <<  |n  |,  (iii)  charge  and  current  neutral¬ 
ity,  (iv)  t)  <<  1  (n  represents  the  fraction  of  incoming  ions 
which  are  reflected) .  This  results  in  the  expression  for  the 
temporal  growth  rate, 
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*2  v ( v  -  kvu  +  ni)2  fti 

Y  =  nQ|k| (V  -  kVu) (V  -  kVu  +  5ni) 


kv. 

J  dVj_  Vj.2 

o 

[(X  ’  — )  3V[  + 

j.  d 

v  Tv~ 

z  . 

fr(vz'vj.,z) 

(1) 


where  Vu  is  the  velocity  of  the  incoming  ion  stream,  f is  the 
distribution  of  the  reflected  ions  limited  to  the  values  of  vz 
such  that  -v  <  vz  <  0#  and  v  is  the  real  frequency  defined  by 


x2  - 


x  Ma  +  1 
a 


JL _ 

X  M_ 

a 


o 


(2) 


with  x  *  k  C/tl.  and  «=  V „/C_.  As  stated  previously,  it  has 
been  shown  that  y  is  positive  for  M_  >  1.5.  For  lower  Mach 

a 

numbers  the  wave  resonant  velocity,  (v  +  ft)/k,  is  such  that  no 
reflected  particles  resonate.  Please  see  Jackson  and  Golden 
[1979]  for  details. 

For  the  case  where  weak  spacial  growth  occurs,  the  wave 
vector,  k,  is  complex  and  the  wave  frequency  real.  Defining 
k(w)  ■  kT(<a)  +  i  k^(w)  and  assuming  |k^|  <<  |kr|,  the  expression 
for  k^  in  terms  of  the  previously  determined  y  is 
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k^(uj) 


(v,k) 


(v,k) 


g 


V  =  (i) 

k=kr(o») 


and  Vg  is  the  wave  group  velocity.  Note  that  in  order  for  k^ 
to  be  small  v^  must  be  finite.  The  linear  analysis  indicates 
that  unstable  waves  occur  for  waves  with  (to,  kr)  in  the  range, 

0  >  (to  +  n.)/k_  >  -  v  .  For  sufficiently  large  M_  (M=  >  2.6 
for  the  ion  acoustic  soliton  with  trapped  electrons)  the  portion 
of  the  (u>,  k^.)  dispersion  curve  which  contains  zero  group 
velocity  whistlers  drops  out  of  this  range.  Two  bands  occur 
in  the  spectrum  of  unstable  whistlers,  one  having  negative 
group  velocities  and  one  with  positive  group  velocities.  The 
whistlers  in  the  first  are  excited  at  z  =  0”  and  grow  while 
propagating  in  the  -z  direction.  The  second  band  of  whistlers 
are  excited  upstream  and  grow  while  being  convected  downstream 
through  the  reflected  beam  to  the  point  z  =  0~.  These  waves 
occur  at  smaller  values  of  |kr|  than  the  negative  moving  waves 
and  will  have  larger  growth  rates  (since  the  expression  for  y 
has  a  |kr|  in  the  denominator).  Therefore,  the  upstream  moving 
modes  will  be  neglected  in  the  following  analysis. 

3 .  Second  Order  Effects 

In  this  section,  the  effect  of  growing  whistler  modes  on 
the  background  plasma  is  analyzed.  The  development  is  analogous 
to  the  second  order  theory  of  Gary  and  Feldman  [1978]  except 
that  here  the  spacial  evolution  of  a  time  average  distribution 
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is  examined  whereas  Gary  analyzes  the  temporal  growth  of  a 
spacially  averaged  distribution.  There  are  some  interesting 
differences. 

Starting  from  the  kinetic  equation  in  one  spacial  dimension. 


Tj.  +  v  -L  +  - i  (e  + 

at  z  az  m^  \  c  j 


f j (Z,v,t) 


and  Maxwell's  equations, 


z  x  E 
az 


1  3B 
c  at 


z  ■—  x  B 

az 


■  4,  I  [  a3v  qj  v  f j  +  J  || 


we  assume  that 


fj(z»v,t)  =  <fj(z,v)>  +  6  fj(Z,V,t) 


E(z, t)  *  <E(Z)  >  +  6  E(z,t) 


STST,  >T 


' ;  *“s 
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B  **  Bqz  +  <B(z)>  +  SB(z,t)  (8) 

where  the  brackets  indicate  time  average.  Using  (6) -(8) 
and  time  averaging  ( 1) — ( 3)  we  arrive  at 

<f j (z,v) > 

=  ~!i  «sf  +  .  j.  it  .>  (9) 

<Ez(z)>  *  l  4tt  qj<f  j  (z ,v)  >  (10) 


where,  as  Gary  points  out,  the  second  order  field  term,  <E  > 

z 

must  be  retained.  Note  that  the  perpendicular  components  of 
<E>  are  zero  if  so  upstream.  Furthermore,  since  the  right 
hand  side  of  (9)  turns  out  to  be  $  independent  ($  is  the 
cylindrical  velocity  variable) ,  it  is  also  consistent  to  set 
<B>  equal  to  zero. 
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One  determines  <5fj,  6E,  k(u)  in  the  usual  manner  by 

subtracting  (9)  from  (3)  and  neglecting  terms  like  6f6E, 

<E  >6f,  and  higher  order  terms.  Fourier  transforming  the 
z 

result  with  respect  to  time  results  in 


[-1“  +  vz  h  -sj  5?]  2-  “> 


-  1 
m 

3 


2-  j<5E(z,u> 


)  +  v  x  6B  ( z ,  g}) 


—  <f^ (Z,V) >. 

3v  J 


(11) 


The  turbulent  fields  are  assumed  to  have  the  form. 


£E(z,u) 


+ 

SE  (o' 


,<d) 


rz 


z 

(x  +  iy)exp[  j  ik(u,o))du],  (12) 


+ 

where  5E*"(o~,<i>)  is  the  spectrum  of  whistlers  at  z  *  o”.  The 
reader  is  reminded  that  these  whistlers  are  excited  upstream 
(where  the  plasma  is  assumed  to  be  in  a  marginally  stable 
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state)  and  grow  as  they  are  convected  back  toward  z  -  o  . 
Using  (12)  in  (11),  solving  for  6f  and  combining  the  re- 

j 

suits  with  Maxwell's  equations  gives  the  dispersion  rela¬ 
tion  for  k(z,to) 


+ 

D~(k,co) 


.22  to  . 

js_  c  -  _  y  _E1_ 

2  2  2n  .to 

to  D  D 


kv  a  kv  . 

,3  z 

d  v  vA  - - - 3— - 

u  -  kv  +  a , 
z  3 


=  0 


(13) 


which  is  valid  for  regions  for  enough  upstream  from  o  that 
ballistic  terms  phase  mix  out.  This  requires  that  the  spread 
in  reflected  beam  velocity  be  large  enough  so  that 


|ki’1(u  +  0j)A(^-)  |»| . 


(14) 


The  approximate  solution  of  (13)  for  k(o>,z)  has  been  discussed 
in  the  previous  section. 


R.H.S.  of  (9) 
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to  -  k  v  +  £2 . 
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Rj (v,z) 
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where 


k+  (-to)  =  k"*(<o)  and  i  6E”  (u>)  | 2  =  |6E+(-co)|2  ( 


have  been  used  (+,  -  refers  to  equation  13.),  the  fields  are 
assumed  axissymetric  and 


12 


This  last  inequality  is  the  usual  requirement  on  the  spectral 

width  for  the  use  of  quasilinear  theory.  It  can  be  shown 

that  both  (14)  and  (17)  can  be  satisfied  if  the  incoming 

ion  stream  is  roughly  a  maxwellian  distribution  and 

2 

S^(B^  =  k  8t  ni/E0  ^  suf f iciently  small.  The  turbu¬ 

lence  is  assumed  weak  enough  so  that  particle  trapping  and 
mode  coupling  [Jackson,  1980a]  can  be  neglected. 

Equations  (9) ,  (10) ,  (12) ,  (13)  and  (15)  can  now  be 
used  to  determine  the  evolution  of  <f^>.  The  initial  con¬ 
dition  at  z=o~  on  is  roughly  specified  by  figure  2. 

The  spectrum  6E(z=o  ,w)  must  be  such  that  <f^(z,v)>  becomes 
marginally  stable  at  some  point,  zu,  upstream  where 
6E(zu,uj)  is  the  noise  spectrum  which  excites  the  instability. 
<5E  (o  ,m)  is  not  actually  determined  in  this  work. 

We  consider  the  change  in  <fj>  just  upstream  of  z=o~, 
far  enough  so  that  ballistic  terms  can  be  neglected,  but 
close  enough  so  that  <fj{z,v)>  can  be  approximated  by 
<fj(o~,v)>  for  all  terms  in  (9)  but  the  first.  This  is 
similar  to  the  procedure  used  by  Gary  (1978). 

The  resulting  equations  are: 
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|  5E-(2,u>)  j 2  =  |6E~(o~,(o)  J2  exp[-2ki“(o“  u)  z]  (19) 


D  [k“(o”,u>),  a),  <  fj(o”,  v)  >]  =  0  (20) 


Equations  (21)  and  (23)  result  from  dividing  (18)  by  v  , 

z 

integrating  over  velocity,  combining  the  resulting  equation 

for  the  density  variation  with  (10)  and  assuming  k.  <<  k 

jL  s 

By  taking  velocity  moments,  some  insight  into  the 
effects  of  the  turbulence  on  the  different  particle  streams 
can  be  obtained.  We  make  the  following  assumptions. 
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where  Vj  and  Cj  are  the  average  velocity  and  thermal  spread, 
respectively,  of  the  jth  stream(The  u,  r,  e  subscripts  indi¬ 
cate  the  incoming  ion,  reflected  ion,  and  electron  streams) . 
Also  used  are  the  assumptions  mentioned  previously, 

[(a*  +  e)Ar  “  Vu  e^  >  >  cu,e  and  l“i  <<:  lnel  *  Note  that 

2  ’ 
kg  simplifies  to 


To  get  the  variation  in  density,  divide  (18)  by  v  t  inte- 

z 

grate,  use  the  assumptions  mentioned  above,  and  use  the  dispersion 
relation,  (20)  .  The  resulting  equations  are, 
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and  W(z)  is  considered  an  integral  operator. 

Note,  that  equations  (25)  through  (27)  combine  to  insure 
charge  neutrality.  In  contrast  to  the  spacially  homogeneous 
case,  charge  neutrality  is  not  automatically  satisfied  by 
(18) ?  it  is  a  result  of  the  long  wavelength,  low  frequency 
characteristics  of  the  turbulence.  Particle  flux  is  auto¬ 
matically  conserved  by  (18)  which  leads  to 
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In  order  to  determine  the  implications  of  (25) -(27), 
we  assume  Vg  to  be  roughly  id  independent  and  consider  some 
characteristic  values  which  occur  if  the  ion  acoustic 
soliton  with  electron  trapping  is  used  as  a  reflecting 
potential, 

dW/dz  >  o , 

Ma  »  2.8,  Mr  =  2.7,  Vg  a  Ma  -  1  •  ( 


This  last  value  for  the  group  velocity  is  a  result  of 
examining  the  linear  solutions  for  k(w)  when  Ma  is  greater 

a 

than  2.6.  It  is  equivalent  to  Gary's [197 8]  assumption, 
a)  -  kVu  <<  ft ; ,  however  this  assumption  was  not  used  in 
arriving  at  (25)-(29). 

It  is  now  evident  from  (25) -(29)  that  the  incoming  ion 
stream  is  slowed  by  the  turbulence  as  it  approaches  z=o” 
from  upstream.  The  electron  density  increases  in  order  to 
preserve  quasi-charge  neutrality.  Reflected  ions  are 
slowed  as  they  move  upstream  from  z=o”  . 
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Further  information  can  be  obtained  by  investigating 
the  momentum  flux  change  of  each  j  stream  defined  by 
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After  some  algebra  [Jackson, 19 81a]  we  arrive  at  the  expres 
sions 
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(33) 


where  the  density  gradients  represent  expressions  (25) -(27). 
The  momentum  flux  gradients  should  be  interpreted  as  the 
force  of  the  turbulent  whistler  fields  on  each  stream.  The 
effect  of  the  second  order  field  <EZ>  is  included  as  part 
of  the  medium  through  which  the  turbulent  fields  act  on  these 
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streams.  In  figure  3  the  relative  size  of  these  forces  is 
illustrated  schematically  for  the  values  in  (29)  .  The  vec¬ 
tors  indicate  the  magnitude  and  direction  of  each  species' 
force  on  the  turbulent  fields  (which  is  represented  by  the 
block) .  Due  to  conservation  of  momentum  flux,  the  forces 
combine  with  the  inertial  force  of  the  block  (the  increase 
in  field  momentum  flux)  to  equal  zero.  The  dashed  vector 
indicates  the  relative  increase  in  field  momentum  flux. 

Note  that  the  change  in  both  field  momentum  flux  and  the 
electron  momentum  flux  (pressure)  is  small  relative  to  the 
change  in  each  ion  stream's  momentum  flux.  The  turbulent 
field  therefore  acts  as  a  medium  for  transmission  of  momentum 
between  the  incoming  and  outgoing  ion  streams.  This  is  the 
same  conclusion  which  Gary  [1978]  made  for  the  homogeneous 
time  varying  case.  In  contrast  to  Gary’s  work,  an  electron 
pressure  force  occurs.  This  force  is  small  but  not  as  small 
as  the  small  parameters,  (23) . 

Using  equations  (31)-(33)  and  assumptions  (23),  it  is 

also  possible  to  show  that  the  change  in  the  z  directed 

velocity  spread  (<5T.  )  is  negligible  [Jackson,  1981a]. 

j  * 

Further  information  can  be  obtained  by  taking  the  mo¬ 
ments, 

f  d3v  ^  (vj,  vj  +  vz2)vz  ±  <£j(V,2)>  , 


) 


which  indicate  the  change  in  perpendicular  and  total  energy 
flux  of  the  jth  species.  Simplified  expressions  (analogous 
to  equations  (25) -(27))  have  been  derived  in  [Jackson,  1981a] 
Using  these  expressions  for  the  present  case  where  the  domi¬ 
nant  unstable  modes  have  positive  group  velocities,  the 
following  conclusions  result: 

(i)  The  total  beam  energy  flux  is  made  up  of  two  parts, 
the  directed  beam  energy  flux,  nrituVr3/2,  which  decreases  in 
magnitude  moving  upstream  from  z=o~  and  the  perpendicular 
energy  flux,  n  Vr<T  ,  which  increases  in  magnitude  mov¬ 
ing  upstream  from  z=o”.  The  increase  in  the  latter  more 
than  offsets  the  decrease  in  the  former  to  give  a  net  in¬ 
crease  in  energy  flux  magnitude  for  the  beam  as  it  moves 
from  z=o”  to  a  point  slightly  upstream.  In  other  words,  the 
beam  is  slowed  and  heated  while  moving  upstream  but  gains 
energy  overall.  If,  in  addition,  the  frequency  of  the 
turbulence  is  such  that  u)  <<  fl^,the  net  gain  in  the  total 
beam  energy  flux  is  much  less  than  the  change  in  the  magni¬ 
tude  of  the  two  component  parts. 

(ii)  As  the  incoming  ion  stream  moves  through  the 

A 

turbulence  in  the  positive  z  direction  it's  perpendicular 
"temperature",  T^u,  increases  and  enhances  the  shock  like 
character  of  this  model.  We  note  further  that  the  total 
incoming  ion  energy  flux  decreases. 


(iii)  There  is  no  change  in  the  perpendicular  temper¬ 


ature  of  the  electrons  (T  ) .  However,  there  is  a  change 

in  the  total  electron  energy  flux  which  decreases  between 

z— o—  and  a  point  slightly  upstream  from  z=o”.  Since  there 

is  no  change  in  Tgz  and  T  o  and  since  the  change  in  directed 

energy  is  small  by  a  mass  ratio,  the  total  electron  energy 

is  small  by  a  mass  ratio,  the  total  electron  energy  flux 

must  be  a  result  of  a  change  in  electron  heat  flux.  Using 

the  previously  mentioned  values  for  V_  and  one  can  show 

g  a 

that  the  change  in  ion  energy  flux  is  much  larger  (  25  times) 
than  the  change  in  electron  energy  flux. 

(iv)  A  comparison  between  the  change  in  the  total  energy 
flux  of  the  reflected  beam  (drwr/dz) ,  the  nonresonant  elec¬ 
trons  and  ions  (d(r  +  r  )/dz),  and  the  turbulent  fields 

W  U  W  C 

(drwf/dz) ,  gives  the  following  approximate  ratios, 


drwf 

dz 


d<rwu  +  W 
3z 


dr 

wr 

dz 


Using  the  sample  values  of  and  M0  is  is  evident  that  most 
of  the  decrease  in  total  non  resonant  energy  flux  goes  into 
increasing  the  magnitude  of  the  total  reflected  beam  energy 
flux,  the  remainder  goes  to  increasing  the  field  energy  flux. 
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4.  Conclusion 

The  analysis  presented  in  the  previous  sections  shows 
the  plausibility  of  the  reflected  ion  shock  model.  Using  the 
ion  acoustic  soliton  as  a  reflecting  potential,  the 
analysis  indicates  that  this  model  is  shock-like  if  M 

cl 

is  large  enough  for  downstream-moving  ,j.ow  frequency  whist¬ 
lers  to  occur. 

A  spacially  vary ing , time  independent , 2nd  order  theory 
has  been  formulated  for  the  turbulence  caused  by  the 
counterstreaming  ions  which  are  characteristic  of  this 
model.  The  pertinent  results  indicate  that  the  main 
action  of  the  low  frequency  whistler  (Almost  Alfven) 
turbulence  will  be  to  slow  and  scatter  (in  perpendicular 
directions)  both  the  incoming  and  outgoing  reflected  ions. 

The  change  in  turbulent  field  and  electron  momentum  and 
energy  is  relatively  small  for  large  Mach  numbers.  There 
is  a  net  loss  of  energy  by  the  incoming  ions  which  is 
gained  by  the  outgoing  reflected  ions. 

As  stated  previously  this  second  order  theory  differs 
from  Gary's  in  a  number  of  ways.  For  instance,  the  effect 
of  an  electron  pressure  may  be  properly  estimated  with  this 
type  of  calculation.  Secondly,  the  expressions  for  gradi¬ 
ents  of  density,  momentum  flux  and  energy  flux  are  expressed 
in  terms  of  quantities  such  a  beam  velocity,  main  ion  stream 


velocity  and  wave  group  velocity  which  are  fairly  easy  to 
interpret.  These  expressions  do  not  require  that  turbulent 
frequencies  be  much  less  than  the  ion  cyclotron  frequencies. 

Also,  this  analysis  considers  a  case  where  some  of  the 
ions  (incoming)  see  a  growing  turbulent  field  and  some  of  the 
ions  (reflected)  see  a  declining  turbulent  field.  This  is 
in  contrast  to  the  time  varying  homogeneous  case  where  all 
particles  see  a  growing  field.  In  consequence,  the  inhomo¬ 
geneous  case  shows  a  net  transfer  of  energy  from  main 
stream  ions  to  beam  ions, whereas, in  the  homogeneous  case, 
the  reverse  is  true.  Finally,  the  effect  of  the  2nd  order 
field,  in  this  case,  is  to  cause  the  electrons  to  preserve 
quasi  neutrality. 

Zero  group  velocity  modes  are  eliminated  from  this 
analysis  by  assuming  Mach  numbers  sufficiently  large  ( >2 . 6) . 
If  these  modes  were  excited  they  would  remain  stationary 
and  grow  in  time  until  they  were  large  enough  to  cause  a 
parametric  instability  [Jackson,  1980]  .  The  resulting 
unstable  whistler  and  ion  acoustic  waves  would  move  down¬ 
stream  and  further  enhance  shock  like  heating  and  slowing. 

We  note  that  by  assuming  Mach  numbers  greater  than 
2.6  it  becomes  possible  for  reflected  ions  to  be  turned 
around  by  the  turbulence  and  to  move  downstream  with  the 
incoming  ions.  These  ions  are  not  included  in  this  model 
and  their  effects  should  be  the  subject  of  further  study. 


4.  Application  to  the  Bow  Shock 

The  model  previously  described  has  a  number  of 
characteristics  which  agree  with  observations  of  the 
earth's  bow  shock.  First,  there  is  a  portion  of  the  bow 
shock  which  can  be  roughly  defined  as  parallel*  .  This 
region  is  characterized  by  the  occurrance  of  low  frequency, 
long  wavelength  electromagnetic  turbulence  having  frequen¬ 
cies  less  than  or  equal  to  the  ion  cyclotron  frequency. 

This  agrees  with  the  reflected  ion  model  if  the  small  k 
instability  is  dominant.  Secondly,  backstreaming  ions  have 
been  observed  in  the  solar  wind  upstream  from  the  bow  shock. 
These  ions  form  two  types  of  distributions,  one  which  is 
sharply  peaked  in  energy  and  direction  and  one  (which  occurs 
most  often)  which  has  a  wide  spread  in  energy  and  direction. 
The  occurrence  of  this  last  distribution  and  of  a  slowing 
solar  wind  is  correlated  with  the  observation  of  low 
frequency  turbulence  occurring  in  the  same  area.  This  agrees 
with  the  slowing  and  perpendicular  heating  described  in  the 
previous  section.  Also,  the  parameters  associated  with  the 
solar  wind  in  front  of  the  shock  are  favorable  to  this  model 
Specifically;  6—1,  T_/T.  >  1,  and  M_—  5.  The  solar  wind 
Mach  number  is  somewhat  higher  them  the  3.1  upper  limit  for 
ion  acoustic  soliton  solutions.  However,  this  3.1  limit  may 
be  modified  by  the  occurrance  of  a  number  of  phenomena  such 
as  Landan  damping,  low  frequency  turbulence  or  a  local  slow¬ 
ing  of  the  solar  wind. 


Since  the  bow  shock  is  rarely  an  exactly  parallel  shock, 
it  is  expected  that  the  ion  acoustic  soliton  (used  as  a 
model  ion  reflector)  is  modified  by  coupling  to  electro¬ 
magnetic  modes.  In  previous  work  [Jackson,  1981b],  it  was 
found,  that,  in  a  low  beta  fluid  plasma,  non  parallel 
propagation  results  in  the  ion  acoustic  soliton  becoming 
a  short  wavelength  ion  acoustic  wavetrain  superimposed  on 
a  long  wavelength  slow  acoustic  wave.  For  larger  beta,  it 
is  expected  that  this  long  wavelength  mode  will  become 
partially  electromagnetic. 

Recently  published  analysis  of  the  two  spacecraft 
I.S.E.E.  data  (June,  1981  issue  J.  Geophys.  Res.) 
has  made  a  large  contribution  toward  describing  the  structure 
of  the  foreshock  region  just  upstream  from  the  quasiparallel 
bow  shock.  A  few  characteristics  described  by  these  papers 
are  not  a  part  of  the  simplified  model  discussed  in  this 
paper . 

First  of  all,  the  bow  shock  is  three  dimensional  as 
opposed  to  the  one  dimensional  model  presented  here.  As  a 
result  of  this,  some  of  the  ions  which  are  in  front  oi  the 
quasiparallel  bow  shock  may  have  originally  been  reflected 
from  a  quasi  perpendicular  region  and  E  x  B  drifted  to  the 
pre-parallel  region  [Bonifazi  and  Moreno,  1981a, b] .  Ions 
reflected  from  a  quasi  perpendicular  region  can  be  expected 


to  have  larger  energy  than  those  which  are  reflected  from 
the  quasi-parallel  region  [sonnerup,  1969],  [Paschmann 
et.  al.,  1980].  In  addition  ions  reflected  at  the  quasi 
perpendicular  shock  make  up  a  narrow  layer  on  the  most 
sunward  region  of  the  foreshock  [Bonifazi  and  Moreno, 

1981a, b]  and  form  the  distribution  which  is  classified 
as  "reflected"  (narrow,  peaked  distributions  in  energy  and 
direction) .  This  is  in  contrast  to  the  parallel  shock  model 
presented  here  where  one  would  expect  the  ions  which  are 
farthest  upstream  to  be  very  "diffuse"  (spread  in  energy 
and  direction) .  This  difference  can  be  resolved  by  remember 
ing  that  the  outer  edge  of  the  foreshock  region  toughes  the 
quasi  perpendicular  bow  shock  [Eastman  et.  al.,  1981], 
Therefore,  ions  reflected  at  this  point  may  run  along  the 
foreshock  boundary  -  xhout  being  scattered  by  the  highly 
turbulent  region  in  front  of  the  parallel  shock. 

Next,  we  note  that  the  model  presented  here  allows  only 
very  weak  magnetic  turbulence  J  6B/BQ J  <<  1  whereas  the  low 
frequency  turbulence  associated  with  the  bow  shock  is  quite 
strong  (|6B/BQ|  ~  1)  [Hopp  et.  al.,  1981].  Very  nonlinear 
waveforms  have  been  observed.  This  is  probably  the  result 
of  allowing  more  reflected  particles  (due  to  larger  e^)  and 
greater  relative  streaming  velocity  (due  to  larger  Mach 
numbers  and  quasi-perpendicular  reflection) .  In  addition, 


low  frequency  compressive  density  variations  which  are 
correlated  with  the  low  frequency  electromagnetic  turbulence 
have  been  detected.  This  may  be  due  to  nonlinear  coupling 
[Jackson,  1980]  or  to  slightly  oblique  propagation  angles 
[Gary  et.  al.,  1981].  One  result  of  these  strong  compressive 
effects  may  be  to  explain  the  large  spread  in  parallel 
velocities  which  is  observed  in  the  distribution  of  "diffuse" 
ions,  but  not  predicted  by  the  model  presented  in  this  paper. 
This  large  spread  in  parallel  velocities  permits  some  ions 
to  be  moving  upstream  without  being  reflected  from  the 
electrostatic  barrier.  It  is  resonant  interactions  with 
these  ions  which  cause  excitation  of  (in  addition  to  right 
hand  waves)  waves  which  are  left  hand  polarized  in  the  plasma 
frame  [Sentmann  et.  al.,  1981]. 

Finally,  observations  by  Bonifazi  and  Moreno  [1981a, b] 
indicate  that  diffuse  i^ns  are  slowed  as  they  move  upstream 
which  is  in  agreement  with  our  model.  Also,  they  found  that 
solar  wind  ions  are  slowed  and  heated  by  tnis  turbulence 
with  a  net  loss  of  total  energy.  This  energy  loss,  minus 
a  smaller  increase  in  field  energy  density,  is  of  the  same 
order  as  the  difference  in  energy  between  the  "diffuse"  and 
the  "reflected"  ion  energy  densities.  (The  total  diffuse 
ion  energy  is  larger  than  the  total  "reflected"  ion  energy 
density.)  A  possible  scenario  is  that  ions  which  are  reflect¬ 
ed  from  a  quasi-perpendicular  region  may  travel  out  in  front 
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{sun side)  of  a  region  of  "diffuse"  ions  loosing  momentum 
and  being  heated.  Gradually  they  slow  enough  to  drift 
into  the  diffuse  region  in  front  of  a  more  parallel  shock 
and  may  then  constitute  a  large  portion  of  the  so  called 
"diffuse"  ions.  The  thermal  energy  gained  by  these  ions 
is  at  the  expense  of  the  solar  wind  kinetic  energy.  If 
this  is  the  case,  the  second  order  analysis  developed  in 
this  paper  is  most  appropriate  near  the  outer  region  of 
the  foreshock  where  the  ions  are  of  the  "reflected"  or 
"intermediate"  classes  and  where  the  low  frequency  waves 
have  not  become  extremely  large  and  nonlinear. 

In  conclusion,  our  reflected  ion  shock  model  seems  to 
agree  qualitatively  with  a  number  of  bow  shock  observations. 
However,  the  three  dimensional  characteristics  of  the  bow 
shock  and  its  large  shock  strength  add  effects  which  must 
be  the  subject  of  future  research. 
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Figure  1.  Schematic  of  reflecting  ion  shock  model. 


Figure  2.  Ion  distribution  at  z  =  o  . 

Figure  3.  Schematic  for  momentum  flux  balance.  Dashed 
line  indicates  change  in  field  momentum  flux. 
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